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ABSTRACT 



Analytical studies have been made of laminar film condensation on a horizontal 
elliptical cylinder in a pure saturated vapor under conditions of free and forced convection. 
Estimation of interfacial shear stress was made in two ways: the first involving an asymptotic 
value of the shear stress under conditions of infinite condensation rate and the second based 
on simultaneously solving the two-phase vapor boundary layer and condensate equations. The 
latter approach enables the determination of the vapor boundary layer separation point. For 
the assumption of asymptotic shear stress, effects of surface tension and pressure gradient in 
the condensate film have also been included. At the extremes of eccentricity, corresponding 
to a circular tube and a vertical flat plate, the results are compared with theoretical and 
experimental work of others. Improvement in the condensation heat transfer coefficient was 
found for elliptical tubes under both free and forced convection conditions when compared 
to circular tubes of the same surface area. In the latter case, this improvement was due mainly 
to the reduced drag of the elliptical tube providing a higher vapor velocity for the same 
pressure drop as that across a circular tube. 
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I. INTRODUCTION 



A. MOTIVATION FOR INVESTIGATION OF ELLIPTICAL TUBES 

For many years, condenser design in the United States has been based on the Heat 
Exchange Institute (H.E. I.) or Tubular Exchange Manufacturers Association (T.E.M.A.) 
standards [1,2] which rely on tube-side empirical relations and averaged data to predict 
heat transfer parameters. Since there is a large degree of uncertainty in the accuracy 
of the prediction, the condensers are over designed to ensure sufficient margin of 
reliability. The result is an excessively large, expensive condenser for the desired 
thermal duty rating requirements. The majority of condensers currently used in the 
Navy were designed under these standards. 

Within the last two decades, computer modeling has been used to more accurately 
predict the heat-transfer coefficients of condensers. There are two basic approaches to 
computer modelling, an integrated and a differential approach, of which the former is 
most commonly used in condenser design (Walker [3]). Integrated approaches divide the 
condenser into zones with associated mean heat transfer and fluid properties which are 
typically determined from semi-empirical relations. These zones are integrated over the 
whole condenser and provide overall performance predictions with respect to shell side 
pressure drop and duty rating. Though easy to use, integrated methods are restrictive 
when used in conditions which are outside of the fluid parameter and geometry 
constraints of the empirical relations employed in the approach. Differential 

approaches solve the basic equations of fluid mechanics and heat transfer and provide 
an understanding of local variations within the condenser, but are too complex for use 
in condenser design. The improvements in accuracy as a rsult of computer design 
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methods are a function of the empirical or theoretical models used in their development. 
Marto [4] provides a comprehensive overview of the state of computer design methods. 

The purpose of empirical and theoretical research into condenser performance is 
to more accurately predict condenser characteristics such that the designer can develop 
reliable, smaller, and hence less costly condensers which meet the same thermal duty 
requirements. Analyses of single tubes allow isolation of individual factors which 
affect condensation without the complexities which result from vapor flow over tube 
bundles or condensate inundation. Since the pioneering work of Nusselt [5], a 
significant amount of theoretical work has been done analyzing laminar film 
condensation on horizontal circular tubes (see Rose [6]), the results of which have 
greatly contributed to the ability to reasonably predict single tube heat transfer 
performance. 

Over the last two decades, major effort has been expended to study effects which 
enhance laminar film condensation heat transfer over that obtained from a plain 
circular tube. The majority of these techniques have focussed on controlling the 
thickness of the condensate film as this is the major resistance to heat transfer. The 
condensate film thickness generally increases with streamwise distance from the top of 
the tube and is dependent on both the rate of condensation and the interfacial shear 
between the vapor boundary-layer and condensate film. Additionally, at high vapor 
velocities, vapor boundary-layer separation occurs which results in a rapid thickening 
of the condensate film downstream of the separation point. Enhancement techniques 
include, but are not limited to, extended surfaces (fins), profiled tube surfaces (roped 
or corrugated tubes), and non-circular tube geometries. Marto [7], Bergles et al. [8] and 
Webb [9] provide comprehensive reviews on such enhancement techniques. These 
techniques enhance heat transfer through an increase in surface area to volume ratio 
and/or the use of surface tension to thin the condensate film. 
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Non-circular geometries contribute to a thinning of the condensate film by placing 
more of the surface in line with the direction of gravity; thus, making the average 
streamwise gravity component in the momentum balance of the condensate film larger 
than that obtained from a circular geometry. Dhir and Licnhard [10] applied a simple 
Nusselt type analysis to an arbitrary plane or axisymmetric body in which the 
streamwise gravity component varied with streamwise length. The results of this 
analysis gave expressions for condensate film thickness and Nusselt number which were 
identical to those derived by Nusselt except that the gravitational constant, g , is 
replaced by an effective constant, g eff , which is defined by 

x (g K ) 413 






ryi'R^dx 

Jo 



(M) 



where x is the streamwise distance from the leading edge of the surface and R is the 
radius of curvature of the axisymmetric body. These expressions can only be applied 
to systems for which curvilinear coordinates are applicable, i.e., for bodies in which the 
radii of curvature are much greater than the film thickness. 

Shklover and co-workers [11,12] analyzed a horizontal cylinder with a logarithmic 
spiral type of surface curvature (see Figure I-l)such that more of the body surface is 
aligned in the direction of gravity and the radius of curvature is continuously 
increasing along the streamwise direction. They manufactured this tube by "mechanical 
deformation of a circular tube." The effect of this profile was to increase the 
streamwise effective gravity component over that found for a circular cylinders and 
to utilize surface tension effects on the pressure gradient to increase the mean film 
velocity, u(x,y), in the streamwise direction. Their analysis considered viscous, gravity 
and surface tension forces in the momentum balance of the condensate film and 
neglected vapor shear at the condensate/vapor interface. These effects resulted in a 
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Figure 1-1. Logarithmic Spiral Tube from Reference [12]. 



thinning of the condensate film and an increase of 20-30% in the overall heat- transfer 
coefficient when compared to a circular tube of the same surface area. These results 
agreed with experimental data. In analyzing the individual effects of gravity and 
surface tension, it was determined that surface tension was most significant over the 
top portion of the tube where dR/dx is largest and accounted for 10-15% of this 
increase. 

An elliptical tube with major axis aligned with the direction of gravity has a 
larger effective gravity than a circular cylinder, as well as an increasing radius of 
curvature over the top half of the tube. In addition, they are more practical from a 
manufacturing standpoint than the non-circular tubes analyzed above. Vapor 
boundary-layer separation should also occur at a point further downstream as compared 
to a circular tube and thus the rapid thickening of the condensate film would be 
delayed. Methods of manufacturing elliptical tubes, and associated problems will be 
discussed in Chapters IV and VI. 
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Wallis [13] conducted flow visualizations (Figure I-2)of water around circular and 



elliptical tube bundles. These photographs provide a qualitative feel for the flow 




Figure 1-2. Flow Visualization, Water over Circular and Elliptical Tube Banks, 
from Reference [13]. 



characteristics around each type of tube bundle. In examining the photographs, it can 
be concluded that the separation point for the elliptical tube bank shifts downstream 
compared to the circular tube bank and results in smaller wake regions behind 
individual tubes. These effects should contribute to better heat transfer performance 
for the elliptical tube bank. Experimental studies by Joyner and Palmer [14] show that 
single phase flow resistance and pressure drop are significantly smaller for elliptical 
tube bundles as compared to circular tube bundles. In power condensers, shell-side 
(vapor) pressure drop has an effect on the overall efficiency of the power plant and on 
the thermal driving force for condensation (T sat -T wa „). Minimizing the vapor pressure 
drop improves power plant efficiency and maintains a relatively constant T sat (for a 





large pressure drop, T $at is reduced as vapor flows through the condenser resulting in 
less thermal driving force). 

Studies involving single-phase heat exchangers and condenser-evaporators indicate 
that elliptical tube heat transfer performance is superior to a circular tube of 
comparable surface area. Ota and Nishiyama [15] conducted an experimental investi- 
gation of single phase forced convection heat transfer of air over an elliptical cylinder 
of minor-to-major axis ratio of 1:3 at various angles of attack. They concluded that 
elliptical cylinders gave improved heat transfer performance over circular cylinders at 
all angles of attack (as a result of fluid turbulence) in the range of Reynolds numbers 
studied (8000 - 79000). Moalem and Sideman [16] conducted a theoretical analysis of a 
horizontal elliptical tube in a condenser-evaporator used in desalinization plants. Their 
study showed a 10 -20% enhancement for elliptical tubes as compared to circular tubes 
with the maximum enhancement achieved at a minor-to-major axis ratio (aspect ratio) 
of 1:4 where the major axis is aligned with the vertical. Huang and Mayinger [17] 
conducted an experimental free convection heat transfer study around elliptical tubes 
and found optimum improvement in heat exchanger performance for tubes with major 
axis aligned vertically and with an aspect ratio of 1:2. Merker and Bahr [18] used an 
analogy between momentum and heat transport to derive a semi-empirical relation from 
which to determine mass transfer rate from the fluid boundary- layer into the free 
stream. The results of their study showed improved heat transfer performance for heat 
exchangers in which the elliptical tubes were spaced wider in the longitudinal direction 
and more compact in the transverse direction. In all these studies it appears that an 
elliptical geometry improves the heat transfer performance of a heat exchanger. 
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B. SIGNIFICANT THEORETICAL STUDIES ON HORIZONTAL CIRCULAR 



TUBES 

The purpose of this survey is to provide a historical perspective to the 
understanding of laminar film condensation on a horizontal circular cylinder. Only 
those works which have made significant advances in this understanding are considered. 
The methods of some of these analyses will later be applied to formulate and solve the 
governing equations for the case of the horizontal elliptical cylinder. Figure 1-3 
provides the geometrical layout and coordinate systems used in these analyses. 




Nusselt [5] used a simple momentum and energy balance to determine the heat 
transfer properties for condensation on flat plates and horizontal circular cylinders. 
In simplifying these physical laws for solution, he made the following assumptions: 



(1) The tube wall temperature is constant. 

(2) The thermophysical properties of the fluid are constant. 

(3) The temperature at the film/vapor interface is T sat . 

(4) The condensate film thickness is small compared to the radius of the tube. 

(5) Heat transfer in the condensate film is one dimensional in the radial direction 
which implies a linear temperature distribution when the film is very thin. 
Convection in the condensate film is neglected. 

(6) The only forces acting on the condensate element are due to viscosity and 
gravity. Inertial forces are neglected. 
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(7) Flow in the condensate film is laminar. 

(8) The vapor is quiescent with no interfacial shear. 



The governing equations 



for this system become 
d 2 u 

q,— + p, g sin<t> = 0 
dx L 



where m is the condensate mass flux rate across the interface, 
is given by 



( 1 - 2 ) 



The mean heat transfer 



5* = 0.943 



A 8 hft p, sin<t> 
ri f L AT 



1/4 



(I -3a) 



for a flat plate at angle <f> to the horizontal (not valid at <|> = 0) and 



= 0.728 



,3 L 2 

h 8 h fg p| 



1/4 



T1 Z D AT 



(I -3b) 



for a horizontal circular cylinder. 

Sparrow and Gregg [19,20] extended the Nusselt [5] analysis by applying boundary 
layer theory to the condensate film. This study incorporated convection and inertia in 
the energy and momentum balance while continuing to neglect interfacial shear. They 
determined that for practical engineering fluids, inertia and convection had negligible 
effect. 

Shekriladze and Gomelauri [21] included interfacial shear by considering 
momentum transfer across the interface due to the condensation process. They used an 
asymptotic expression for the interfacial shear based on an infinite condensation rate 
given by 

= MU* - «*) ( N4 ) 

where is the streamwise velocity of the vapor at the outer edge of the boundary- 
layer. It was also assumed that U^»u 6 such that the film velocity may be neglected and 
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that inertia and convection effects in the condensate film could also be neglected. As 
noted by Rose [6], the validity of neglecting u 5 requires that A^AT/r^h^ be small (large 
Prandtl number). For forced convection heat transfer on a horizontal circular tube, 
Shekriladze and Gomelauri determined the following expression: 

Nu = 0.9 t 1 ' 5 ) 



where 



Re 



P , U m D 
n, 



( 1 - 6 ) 



is the two-phase Reynolds number. When the effect of gravity cannot be neglected 
(mixed convection) the mean heat-transfer coefficient is given by 

Nu = OMs/tiefl +V1 + 169 F (, ' 7) 

where 

F = n/ ^ - - (1-8) 

X l AT Ul 

and measures the relative effects of vapor velocity and gravity. It should be noted that 
based on potential flow outside the vapor boundary layer, the interfacial shear will 
always be positive and therefore this method does not predict vapor separation as would 
occur in reality. Based on a minimum separation angle of 82° for single phase flow over 
a circular cylinder without suction and no heat transfer after separation, Shekriladze 
and Gomelauri conservatively estimated that the heat-transfer coefficient would be 
reduced by about 35%. Actual separation occurs at a point further downstream and 
results in a heat-transfer coefficient between the result for separation at 82° and no 
separation. These expressions have been found to be reasonably accurate for cases of 
high condensation rate and are useful for their simplicity. 

Fujii et al. [22] conducted an analysis of mixed convection condensation on 
horizontal cylinders using the two-phase boundary -layer equations. Their formulation 
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neglected inertia, convection and pressure gradient effects in the condensate film. The 
governing equations for the model are given by 



dU + dV 
dx dy 



0 



U— + V— = U 
dx dy * dx p v dy 2 



(I -9a) 



for the vapor boundary-layer and 

du dv 



+ — - = 0 
dx dy 



n #u . . 

0 = t), — - + p, g sm(J) 



(I -9b) 



m 



= P^f'udy = 

dx J o 



X,AT 



6 h 



fg 



for the condensate film. The compatibility equations at the interface are given by 



U t =0 



i du 
1,1 * 



= •n, 



/>=6 



in 



= t. 



Pi 



( as 
1 ax 



/ y=6 
= -Pv V y-i 



(I-9C) 






An approximate integral solution of the momentum equation for flow over a cylinder, 
with suction by Truckenbrodt [23], was modified to more closely agree with the 
numerical solutions of Terril [24]. This method enabled the determination of the 
interfacial shear and subsequent solution of the governing equations for the condensate 
film. This technique also predicted the point of vapor separation and subsequently, a 
more accurate heat- transfer coefficient. For locations downstream of the vapor 
boundary-layer separation point, the interfacial shear was assumed to be negligible and 
simple Nusselt theory used to determine the heat transfer. Local and mean Nusscit 
numbers were numerically determined for boundary conditions of uniform wall 
temperature and uniform wall heat flux. Their results were in fair agreement with 
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with experimental data. In addition to the dimensionless parameter, F, a new 
dimensionless parameter, G, was defined as 

G = — k ‘ P< 'Hi (I- 10) 

h/gi p v n v 

which is a measure of the condensation rate. As G-+<», the results of Fujii et al. [22] 
agreed with the previous results of Shekriladze and Gomelauri [21]. Lee and Rose [25] 
noted an error in the formulation of Fujii et al. [22] in that the interfacial shear terms 
in the condensate film momentum equation should have been divided by 2. They re- 
analyzed the system of equations using both the modified and unmodified method of 
Truckenbrodt [23] and found little difference whether one used the modified or 
unmodified method of Truckenbrodt. 

To date, Gaddis [26] has conducted the most comprehensive study of condensation 
on horizontal circular tubes using the two-phase boundary-layer equations. His analysis 
neglected surface tension in the momentum equation as well as viscous dissipation and 
pressure in the energy equation of the condensate film. He determined that inertia 
effects were negligible for most media (with the exception of liquid metals) and that 
convection effects were somewhat significant for viscous liquids with Pr :» 1, similar 
to the conclusions of the simpler analyses of Sparrow and Gregg [19, 20]. He identified 
three regimes of behavior for flowing vapors. For low vapor Reynolds numbers, Re v , 
gravity effects dominate the heat transfer and Nusselt’s analysis adequately represents 
the heat transfer behavior. For high Re w , condensation was shear controlled and the 
Nusselt number is proportional to the square root of Re v as found by Shekriladze and 
Gomelauri [21]. The intermediate region, in which gravity and vapor shear are both 
significant, results in no simple relation. Gaddis also analyzed several cases of flow 
separation (vapor boundary-layer and/or condensate film separation). For low Re w ( * 
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10 2 ), vapor boundary-layer separation occurred with no flow reversal in the liquid film. 
This result is expected as the pressure gradient effect on the condensate film, which 
tends to slow the condensate flow over the rear half of the tube, is a function of the 
vapor potential velocity. At moderate Re v (~ 10 3 ) vapor boundary-layer separation nd 
condensate flow reversal near the interface occurred with a resulting decrease in the 
mean fluid velocity and rapid thickening of the condensate film. At high Re w (~ 10 5 ), 
no vapor boundary- layer separation occurred but the condensate film separated at the 
wall causing a sharp reduction in mean condensate velocity with a rapid thickening of 
the film as a result of a significant, adverse pressure gradient. Fluid separation at their 
respective boundaries is defined as the condition at which the shear is less than or equal 
to zero. 

diMarzo and Casarella [27] conducted a similar analysis using a more general 
solution technique and arrived at the same results as Gaddis [26]. In their analysis of 
the flow separation phenomenon, they provided guidance on determining heat transfer 
performance once separation had occurred. At low Re v , gravity effects dominate and 
no flow reversals or thickening of the condensate film occurred. The use of a Nusselt 
type analysis with no interfacial shear would seem prudent in the region beyond vapor 
separation. At moderate to large Re w , pressure gradient effects dominate. The adverse 
pressure gradient over the back half of the tube causes flow reversal and a rapid 
thickening of the condensate film. In this case, it would be prudent to neglect heat 
transfer completely beyond the separation point of the vapor-boundary layer or 
condensate film. 

Rose [28] studied the effects of pressure gradient in the condensate film. In this 
case, the pressure gradient on the condensate film is due to the pressure gradient of the 
vapor, as determined by potential theory, which is impressed on the condensate film. 
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He simplified the analysis by using the Shekriladze- Gomelauri model [21] which ignores 



inertia and convection in the condensate film and uses the asymptotic expression for 
the interfacial shear. As a result of this analysis, a dimensionless parameter, P , was 
defined by 



which represents the pressure gradient effect (Gaddis [26] had a similar combination 
of dimensionless parameters which were equivalent to P). When P = 0, the governing 
equations reduce to those of Shekriladze and Gomelauri [21]. He concluded that the 
effect of including pressure gradient was to improve heat transfer over the forward 
half of the tube since the pressure gradient over this region is favorable (tends to 
increase the mean film velocity which reduces the condensate film thickness). As a 
result of the formulation, he found that for cases where P > F/8, a critical angle at 
some point on the rear of the tube was reached where d6 /d<J) - «. In this case, it was 
not possible to obtain a solution over the entire tube. It was postulated that this critical 
angle might indicate some instability followed by some degree of waviness or 
turbulence. It was noted that this critical point was reached prior to the point at which 
the condensate film separated. For conditions which permitted solution over the entire 
tube (i.e. P < F/8), any increase in heat transfer achieved over the forward part of the 
tube was balanced by a decrease in heat transfer as a result of the adverse pressure 
gradient over the back half of the tube such that there was little change in the mean 
heat- transfer coefficient. Since the pressure gradient due to potential flow of vapor 
just outside the vapor boundary-layer is given by 



P = Py h h 
P, at 



( 1 - 11 ) 




( 1 - 12 ) 
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which is symmetric around the surface of the tube from front to back, the favorable 
effect over the front half of the tube is exactly offset by the adverse effect over the 
back half of the tube. Referring to the studies conducted by Gaddis [26] and diMarzo- 
Casarella [27], in the region of moderate Re y , condensate film separation from the tube 
wall does not occur but the pressure gradient sufficiently retards the film such that 
flow reversals occur and the film still rapidly thickens. In the present study, the effects 
of pressure gradient using the analysis of Rose [28] resulted in a condensate velocity 
distribution which approached zero as the critical angle was reached. This phenomenon 
may be a result of the model used in formulating the problem. Since Shekriladze and 
Gomelauri [21] used an asymptotic expression for the interfacial shear, negative 
velocities cannot exist in the velocity profile. The positive velocity distribution is 
retarded by the adverse pressure gradient. In the case of P> F/8, the pressure gradient 
is large enough to cause complete stoppage of film flow at <f) c and the film thickness 
increases rapidly. 

Krupiczka [29] examined the effects of surface tension due to film curvature on 
condensation on circular cylinders. He used a simple Nusselt type model but included 
surface tension in the momentum equation. In his development he did not assume the 
film thickness was much less than the radius of the cylinder until after the inclusion 
of the surface tension term to account for the curvature of the film. The resulting 
equation was a second order ordinary differential equation which required two initial 
conditions. The first condition was given by the symmetry of the problem. However, 
the initial thickness was not obvious and was arbitrarily chosen to be that obtained 
from the Nusselt model. He concluded that the effect of surface tension was small on 
the forward part of the tube and increases in significance over the back part of the tube 
due to the rapidly changing film thickness. This significance was dependent on the 
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magnitude of dimensionless parameter, A, given by 

A = -We 
4 

which is dependent on the surface tension of the fluid through a "modified ,, Weber 
number (We = a/(pgr“)) and the radius of the cylinder, r . The conclusion was that 
surface tension can be significant for small diameter tubes and wires where, as r 
becomes smaller, the film thickness becomes relatively significant and hence the film 
curvature cannot be neglected. For most practical cases, and therefore the surface 
tension effect due to film curvature can be neglected. The effect of assuming the 
initial film thickness to be that given by Nusselt [5] was checked and verified to have 
minimal effect on the mean heat-transfer coefficient. 

C. SIGNIFICANT THEORETICAL STUDIES ON HORIZONTAL ELLIPTICAL 
TUBES 

Some theoretical work has been done on horizontal elliptical cylinders in a 
quiescent pure vapor using a Nusselt type model. The physical orientation and 
coordinate system is shown in Figure 1-4. 

Cheng and Tao [30] approximated the surface of an ellipse by several circular arcs. 
They analyzed condensation on these arcs using the same assumptions as Nusselt. The 
ellipse was aligned such that the major axis, a, was aligned with the direction of 
gravity. From their numerical results they determined that the heat- transfer 
coefficient decreased with increasing eccentricity, k (defined as the ratio of the minor 
to major axis). Values of k were varied from that of a vertical flat plate (k = 0) to that 
for a horizontal circular cylinder (k ~ 1). The mean Nusselt number was determined by 
a surface area weighted average of the mean Nusselt number for each circular arc and 



3t| ]ja 

~2 3 D 

Pi 8 r Pr 



1/4 



(M3) 



15 







Figure 1-4. Geometry for Film Condensation on a Horizontal Elliptical Cylinder. 



Cj is an integral constant which links each arc and i is the number of arcs used in the 
approximation. In the practical range of eccentricity (0.3 - 0.6), the mean heat transfer 
coefficient was increased by 10 to 18% over that of a circular tube with the same 
surface area. Ali and McDonald [31] conducted a similar type analysis as Cheng and 
Tao without the circular arc approximation as a first estimate for condensation on 
inclined circular tubes. 

Wang et al. [32] used the Nusselt assumptions to analyze condensation on 
horizontal elliptical tube for which the major axis is oriented at an angle, a, with 
respect to the vertical axis. They obtained an expression for the mean heat- transfer 



is given by 




(1-14) 



where 




( 1 * 15 ) 
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coefficient given by 



where 



a 



2 . 3 

1 h fg Pi 8 h 1 

4 ^ fib 



1/4 




(1-16) 



Z(p) = sin'^p fj 



ab 



fl 2 sin 2 (p-a) + 



fc 2 sin 2 (p-a) 



3/2 

sin 1/3 p # 



(1-17) 



The results of this study showed that the maximum mean heat- transfer coefficient is 
obtained when the major axis is aligned with the vertical (a = 0° ). These theoretical 
results were validated with experimental data for an elliptical tube with semi-minor 
axis of 8 mm and semi-major axis of 22 mm. Sheng and Cha’o [33] noted that the mean 
heat- transfer coefficient was incorrectly determined since the heat- transfer coefficient 
is averaged over the surface area and the radial distance, /*, is not constant for an 
ellipse. In the course of this current study, the problem formulation of Wang et al. [32] 
was rerun with the correct expression for the mean heat- transfer coefficient (as later 
derived in Chapter II). The corrected mean heat- transfer coefficient of the elliptical 
tube analyzed by Wang et al. is 11.3% larger then a circular tube with the same surface 
area. Additionally, it should be noted that the heat transfer performance of this 
elliptical cylinder was better than the circular cylinder for angular orientations up to 
a » 50°. The corrected theoretical results are approximately 5% lower than those 
determined by Wang et al. and more closely agree with their experimental data. Figure 
(1-5) shows the heat transfer enhancement for elliptical tubes of varying eccentricity 
and orientation angle. 

Sheng and Cha’o [33] studied the effects of surface tension and variable wall 
temperature on condensation on a horizontal elliptical tube. The remaining assumptions 
of the analysis were the same as those above. For the wall temperature, a cosine 
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Figure (1-3). Elliptical Tube Mean Nu Compared to Circular Tube Mean Nu for 
varying Eccentricity and Orientation Angle. 

distribution of the form 

AT = (T^-T^Jd - A COS0) (1-18) 

was used (which has been shown (Memory and Rose [34]) to be in good agreement with 
experimental data) where the value of A depends on the ratio of the outside to inside 
heat-transfer coefficient. They determined that variable wall temperature affected 
local values but not the mean values of the heat- transfer coefficient. For the surface 
tension effect, it was assumed that the film thickness was much smaller than thv r adius 
of curvature of the elliptical surface and thus the surface tension effect was due fely 

to the curvature of the tube wall. This is in contrast to the study of Krupiczka [ or 

a circular tube where the surface tension was due to the curvature of the conde . *i te 
film. Over the back half of the tube the radius of curvature decreases with stream rise 
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distance resulting in a retarding effect on film flow. As the flow is slowed, the film 
thickness increases rapidly and at some critical angle becomes infinite in magnitude. 
This phenomenon is similar to that described by Rose [28] for the pressure gradient 
effect. The surface tension effect was determined to be significant for k < 0.6 . In 
comparison to a circular tube, surface tension causes a favorable pressure gradient over 
the front half of the tube resulting in a greater film velocity, thinner film thickness and 
improved heat transfer. It has an opposite effect on the back half of the tube which 
tends to negates this improvement. The gravity component of the momentum balance 
is the driving force behind the enhancement in elliptical tubes. Consideration of 
surface tension results in a slight decrease in the mean Nusselt number as compared to 
the situation in which the surface tension is neglected. 

A potential advantage of an elliptical tube compared to a circular tube is the 
difference in vapor flow characteristics as a result of a better streamlined shape. 
Panday [35] developed an explicit numerical method for two dimensional film 
condensation and applied it to the case of downward flowing vapor over elliptical 
cylinders. Convection and inertia were included in the condensate film as well as 
surface tension and pressure gradient (as a result of potential flow of vapor outside the 
vapor boundary-layer). The interfacial shear was approximated using the asymptotic 
expression for infinite condensation rate. Several errors were found in the expressions 
for surface tension and pressure gradient as a result of an incorrect analysis of the 
differential streamwise length, dx. The first error is the result of assuming that dx = 
rd(f>. This relationship assumes that the radial distance from the centroid of the ellipse 
is constant over the interval of the parametric angle. This fact is not true (as will be 
shown in Chapter II) and results in an error whose magnitude is dependent on the step 
size used in the numerical procedure. The second error involves the expression for r. 
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If the parametric angle is measured from the vertical axis of the ellipse (as appears to 
be the case in Panday’s expression for vapor velocity potential), then r should be given 
by 

x = b sin<J> 

y = a cos<J> (1-18) 

r ~ \jx 2 +• y 2 = a\Jcos 2 $ + & 2 sin 2 <|> 

Panday used an expression for radial distance given by 

r = ayjsin 2 <J> + £ 2 cos 2 4> (l - 19) 

The magnitude of the resulting error is dependent on the eccentricity of the ellipse. 
Additionally, there appears to be an error in the gravity component of the momentum 
equation. Panday uses a body force given by pgsinO, which is true only for circular 
geometries but is untrue for elliptical geometries. His results will therefore not account 
for the improved performance which are obtained by placing more of the surface in 
line with the direction of gravity. Panday’s conclusions were that the overall heat 
transfer was reduced for elliptical tubes at low velocities (as compared to a similar 
circular tube) due to a rapid thickening of the film at the stagnation point and the 
overall heat transfer is increased for high velocities due to increased interfacial shear. 
These conclusions are opposite to what one may expect. At low velocities, the film 
thickness near the stagnation point for elliptical cylinders is thinner than the equivalent 
circular tube as a result of increased effective gravity. This effect results in improved 
heat transfer. Pressure gradient effects due to potential flow and interfacial shear 
should be negligible since the velocity is small. At high velocities, the streamlined 
geometry of elliptical cylinders results in lower interfacial shear and pressure gradient 
effects as compared to circular cylinders. These effects should result in a thickening 



20 



of the condensate film as compared to circular cylinders and hence reduced heat 
transfer for the same effective diameter and free stream vapor velocity. 

D. OBJECTIVES OF THIS THEORETICAL STUDY 

The results of previous studies indicate that elliptical tubes can be used to 
increase condensation heat transfer in low vapor velocity condensers as compared to 
circular tubes with the same surface area. The present investigation examines the 
effects of vapor shear, pressure gradient and surface tension in laminar film 
condensation on a single horizontal elliptical tube with its major axis aligned with 
gravity and the free stream velocity. Interfacial shear is estimated using both the 
simple assumptions of Shekriladze and Gomelauri [21] and the more complex technique 
of Fujii et al. [22], The latter case calculates the angle at which the vapor boundary- 
layer separates such that the effect of reduced drag on the mean heat- transfer 
coefficient can be evaluated. 
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II. THEORETICAL DEVELOPMENT 



Following the theoretical developments used for the case of laminar film 
condensation on a horizontal circular tube, theoretical models are developed for laminar 
film condensation on a horizontal elliptical tube of eccentricity, k. This development 
starts with a Nusselt [5] type model in a quiescent vapor, and then adds forced and 
mixed convection as in the models of Shedkriladze and Gomelauri [21] and pressure 
gradient as considered by Rose [28]. The pressure gradient takes into account the 
effects of potential flow outside the vapor boundary-layer as well as surface tension. 
Finally, a model is developed which analyzes the vapor boundary-layer and boundary- 
layer separation following Fujii et al. [22]. Pressure gradient and surface tension are 
not considered in this model due to the complexities introduced by the use of the two- 
phase boundary- layer equation. Where possible these elliptical models are checked 
against existing theories for the "limiting" eccentricities of a circular cylinder (£=0), 
vertical flat plate (/:=!_) and horizontal flat plate {k -* «). 

A. FACTORS RELATED TO ELLIPSE GEOMETRY 

Consider an elliptical cylinder whose cross-section is oriented such that the major- 
axis is aligned with the vertical as shown in Figure II- 1. The eccentricity, k , of the 
ellipse is defined by 

it = ^ . (ii- 1) 

a 

where a and b are the semi-major and semi-minor lengths, respectively. Functions 
related to the geometry of the ellipse are initially developed in a cartesian coordinate 
system (x,,y ,) whose origin coincides with the centroid of the ellipse. This is then 
transformed into a cylindrical coordinate system (r,0) where r is the radial distance 



22 



from the ellipse centroid to a point on the ellipse surface and 0 is the angle measured 



from the vertical. The transformation equations are 

x x = r sind 
y x = r cos0 . 



(II -2) 



The elliptical surface can also be defined by a parametric angle, 4>, measured from 
the upper semi- major axis, such that the y t coordinate on a circle of radius a translates 




Figure II- 1. Geometry for Film Condensation on a Horizontal Elliptical 
Cylinder with Major Axis Parallel to Gravity. 



to the y ! coordinate on the ellipse and the x 1 coordinate on a circle of radius b translates 

to the x 1 coordinate on the ellipse. The transformation equations are given by 

x x = b sin<J> 

= a cos<j> . 



(II-3) 



Combining Equations (1-2) and (1-3) results in a relationship between r,0 and 4>: 
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(II -4) 




cos<t> = — cos0 ( b ) 

a 




(c) . 



It is assumed that the condensate thickness, 6, is much smaller than the radius of 
the elliptical surface. Therefore, the equations for the condensate film and vapor 
boundary- layer are developed using an orthogonal, curvilinear coordinate system, (;c,y), 
where x is the streamwise distance along the elliptical surface and y is the distance 
normal to the surface. 

1. Radial Distance, r 



Using the parametric transformation equations (Equation (II - 3)), an equivalent 
expression for the radial distance is given in terms of <t> by: 



In cartesian coordinates, the surface of an ellipse is given by 




(II ' 5 ) 



The radial distance, r, can be determined by 




(II -6) 



r(4>) = a{ cos 2 4> + k 2 sin 2 <}> 



(II -7) 
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2 . 



Radius of Curvature of an Ellipse, R 



The curvature, k^), and the radius of curvature, R(x 1 ), of a function, f(Xj), 



is given by 



= 



[1 + 0 r/ (* 1 )) 2 ] 3/2 



RixJ = 1/kCCj) 



(II -8) 



where the prime denotes differentiation with respect to Xj. Solving for y 1 as a function 
of Xj using Equation (II-5), determining the radius of curvature using Equation (II - 8) 
and transforming the results to parametric coordinates gives the radius of curvature as 
a function of <J>: 

m = - [sin 2 4> + • (n-9) 

k 



3. Streamwise Length, * 

Since the radial distance to a point on the elliptical surface is not constant, 
the streamwise length is not proportional to 0 as it would be for a circle. Additionally, 
as eccentricity varies, the streamwise length for a given 0 is not the same. To obtain an 
expression for x in terms of parametric angle 4b consider first a point on the ellipse 
surface as shown in Figure II-2. Moving a small distance, dx y results in incremental 

changes, dO and dr . The resulting relationship between x,r and 0 is given by: 

(dx? = ( dr ) 2 + (r d&f (II- 10) 



Expressions are now needed for dr and d6 in terms of d<j). 

achieved by differentiating Equation (II-7): 

dr = 1 a (* 2 -D sia2< fr_ d<f> . 

^ v/cos 2 * + £ 2 sin 2 4> 



For dr, this is simply 
(II- 11) 
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Ellipse. 



For d0, taking Equation (II-4a), substituting Equation (II - 7) for r and differentiating 
the resultant expression results in: 

cose de = k - os * . (ii- 12) 

[cos 2 4> + fc 2 sin 2 4>] 3/2 

Combining Equations (II - 4b) and (II - 7), and substituting for cos0 in Equation (11-12) 

results in an expression for dQ in terms of d<f>: 

dQ = — - — d<\> . (II-13) 

cos 2 <}> + k 2 sin 2 4> 

Substituting Equations (II- 11) and (11-13) into Equation (II- 10) gives an expression for 
dx as a function of d<f>: 



dx = 



k 2 +—(k 2 -l) 2 siir2<J) 

a — - — d§ . 

\ cos 2 4> + k 2 sin 2 4> 



(11-14) 



This may be integrated to obtain the streamwise length. For compactness, a function 
X(4>) is defined by 



X(<|>) = a 



\ 



k 2 + — (A 2 - 1)^!!! 2 ^^ 
4 

cos 2 <J> + k 2 sin 2 <J) 



(IM5) 
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in later developments, a characteristic length is used to non-dimensionalize 
the heat transfer parameters in the model. An effective diameter, D e , is defined as the 
diameter of a circular cylinder having the same surface area as an elliptical cylinder 
and will be used as a characteristic length. The effective diameter is given by 

nD = 2 rx(4>)d<t> ■ (II -16) 

* Jo 



Streamwise length is non-dimensionalized by: 

jc* = — 

7 



(11-17) 



This enables direct comparison between elliptical and circular cylinders since it 
represents an equivalent surface area. 

4. Streamwise Gravity Component, g x 

The component of gravity in the streamwise direction is tangent to the 
elliptical surface and is a function of streamwise location. A line tangent to the 
surface, defined by the slope of the surface, is obtained from Equation (II-5), 

fy a * i 



dx, 






(11-18) 



A unit vector in the streamwise direction is then given by 



T = 



1 



N 



a 2 „ 2 



f7~2 2 ? a * 

r *1 1 “ T *1 j 



ib 1 - *i) - ^ 

b l 



(11-19) 



The streamwise gravity component is given by the dot product of the gravity vector and 
the tangent vector, 



8x = 8 




\ 



C b 2 - xt) + 




(11-20) 
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Transforming Equation (11-20) into parametric coordinates results in the expression: 



8 X = g /i(4>) 



where 



W) = 



sin({) 

+ fc 2 cos 2 4> 



(II -2 1 ) 



5. Potential Flow at Ellipse Surface, 

When considering forced convection due to the flow of vapor over the 
condensate film, the velocity of the vapor influences the condensate film thickness 
through the vapor shear at the film/vapor interface. The vapor velocity at the interface 
will be dependent on the vapor velocity outside the vapor boundary-layer. This velocity 
is determined from potential flow theory. Assuming that the film thickness, 6, and the 
vapor boundary-layer thickness, A, are much smaller than the radial distance, r, of the 
elliptical surface (and therefore may be neglected), the potential flow, U^, about an 
elliptical surface with major axis aligned with the vapor free stream velocity, U tt , is 
given by 

= V- ; ■ (l+k) • (II -22) 

yl + k 2 c ot 2 ^ 

Details of the derivation of this expression are provided in Appendix (A). 

6. Ellipse with Major Axis Perpendicular to Gravity 

For the case where the ellipse major axis is perpendicular to the vertical as 
shown in Figure II-3,the previously determined equation for streamwise component of 
gravity remains the same. For 0 < k < 1, the major axis is aligned with gravity. For 1 
</:<«, the major axis is perpendicular to gravity. 
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Figure II-3. Geometry for Condensation on Horizontal Elliptical 
Cylinder with Major Axis Perpendicular to Gravity. 



B. FREE CONVECTION CONDENSATION ON AN ELLIPTICAL CYLINDER 
(NUSSELT [5] TYPE ANALYSIS) 

Consider a quiescent, saturated vapor at a temperature T sat , condensing on a 
horizontal elliptical cylinder of eccentricity, k y and semi- major axis, a, aligned with the 
direction of gravity as in Figure II- 1. The same simplifying assumptions are made as 
proposed by Nusselt in his analysis of condensation on a horizontal circular cylinder. 
These assumptions are: 

(1) The vapor is pure and quiescent. 

(2) The tube wall temperature is uniform and constant. 

(3) The thermophysical properties of the fluid are constant and evaluated at some 
given reference temperature. 

(4) The temperature at the film/vapor interface is T sat . 

(5) The condensate film thickness is small compared to the radius of the ellipse. 

(6) Heat transfer in the condensate film is one dimensional in the radial direction 
providing a linear temperature distribution across the film (based on 6<R). 

(7) The only forces acting on the condensate film are viscous and gravity forces. 

(8) Flow in the liquid film is laminar with no waviness and no vapor boundary 
layer separation. 
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A fluid element of unit width is analyzed in detail and shown in Figure II - 4. 




Conservation of mass for the fluid element requires that, 

d 



K‘“ *] * - *} 



dx + m dx = 0 



where m is the local condensation mass flux rate. 
Defining the mean film velocity as 

u = —f 6 u dy 



(11-23) 



(11-24) 



Equation (11-23) simplifies to 

m = (II ’ 25) 

Conservation of momentum for the film element is a balance of viscous and body 
forces which reduces to 

<1,^-7 *P, </,(♦>■<> - 
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Here, an assumption is made that the fluid is Newtonian, i.e., 



du 

n/ 



(11-27) 



Conservation of energy is a balance of latent heat from condensation and heat 
conduction through the condensate film which reduces to 



V* * = * 



m = 



^7 (^sof Vq//) 



(11-28) 



The boundary conditions for the momentum equation are 

Vo = 0 



t. = 



(du 

'U 



= o 



/y=5 



(11-29) 



Integrating the momentum equation subject to the boundary conditions results in 
a local film velocity of 

vo = ^/,(<t>) [*y - \y 2 ) (h- 30 ) 



and a mean film velocity of 




(11-31) 



Combining the continuity equation, energy equation and expression for film mean 



velocity (Equations (II - 25), (11-28) and (11-31)), results in 



X, AT 



a\ 


\ Pi 8 


dx\ 


[H 



(11-32) 



Using the relationship between dx and dtp from Equations (11-14) and (11-15) and 
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carrying out the differentiation results in 



— - - = Id 3 //( 4>) + 3 /,(<t>) 5 2 — 

h fs 5 3r) f x(<l>) l Jl ' <*t> 



(11-33) 



A dimensionless parameter, Z 1? is defined by 



Z, = — 
1 *1 



(11-34) 



where B , is given by 



= 



X, ti, A T 

pF \ g 



(11-35) 



Note that the equivalent parameter used by Nusselt [3] multiplies Equation (11-35) by 
a factor of three and uses radius in place of effective diameter. Substituting Equations 
(11-34) and (11-35) into (11-33) results in a first order ordinary differential equation for 

z* 

dZ^ + 4 /i / (4>) z = 4 x(4>) (11-36) 

3/j(4)) 1 D /,(4)) 



This equation can be solved exactly by using an integrating factor which would 
require numerical integration or can be solved numerically using a forward stepping or 
propagation technique. The initial condition is determined based on the symmetry of 
the problem which requires that (dZj/d <J>)^ =0 = 0. Solving Equation (11-36) at <t> = 0 
results in the initial condition 



"l $=o 



3 X(0) 

D' /.'(0) 



3 a £ 2 

D < 



(11-37) 
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The non-dimensional film thickness, d , and film thickness, 6, are determined 
respectively by 



6*(4>) = Zj( 4>) 1/4 and 

d(4>) = < Zj(4>) Bj } 1/4 . 

The local heat- transfer coefficient, a, and mean heat- transfer coefficient 
respectively by 



(11-38) 
are given 



a(d>) = — — and 

m 



a 



C a ( *> 

!>s 



-%r /’*(*) X(4>) <*4> , 
tiD •'o 



(11-39) 



and the local Nusselt number and mean Nusselt number by 



Nu (< (>) = 



tt(4>)Q, 



m 



and 




2 fn 1 

11 Jo 5(<{>) 



x(4>) <*► 



(11-40) 



For the case where = 1 (circular tube); f 1 (4>)=sin<t>, f j '(4 > )=cos<t>, x(40 =a and D e = 

2a and Equation (11-36) reduces to 

<&i + 4 cos<t) z = _2_ ' (11-41) 

J4> 3 sin<(> 1 sin<t> 

which is the same as that found by Nusselt [5] except for a difference in the definition 
of B v 

The above development can also be applied to a horizontal elliptical tube whose 
major axis is perpendicular to the direction of gravity as shown in Figure II-3. In this 
case, b is greater than a , i.e. k > 1. 
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C. FORCED CONVECTION CONDENSATION ON AN ELLIPTICAL CYLINDER 
(SHEKRILADZE- GOMELAURI [21] TYPE ANALYSIS) 

Consider a saturated vapor at temperature T sat , flowing downward over a 
horizontal elliptical cylinder with free stream vapor velocity, U ro . The elliptical 
cylinder has eccentricity, and semi-major axis, a , parallel with the direction of 
gravity and vapor flow. The same assumptions as for the case of free convection are 
used here with the exceptions: 

(1) The force of gravity is neglected. 

(3) U^u y=6 . 

(4) The interfacial shear stress is approximated by an asymptotic expression 
assuming an infinite condensation rate, 

= mU + (11-42) 



Conservation of mass and energy are the same as previously derived (Equations 
(11-25) and (11-28)). The only forces acting on the fluid element are viscous forces. 
Thus, conservation of momentum reduces to 

1, - o , (11-43) 

dy i 



with the boundary conditions, 



«,-o = 0 



> 1 / 



BL 



** = m U* . 



(11-44) 



Integrating the momentum equation subject to the boundary conditions results in a 
local velocity of 

m U* /II ICN 

u(y) = * y (H-45) 

’ll 
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and a mean film velocity 



u_ - 



m U* b 

n, 2 



(11-46) 



Using the expression for potential flow over an elliptical cylinder (Equation (II-22))and 
defining a velocity function, f 2 (4>), such that 



M) = 



1 +k 






(11-47) 



the mean film velocity can be expressed as 

m U 



u_ = 



— m\ 
0/ 2 



(11-48) 



Combining the continuity equation, energy equation and expression for mean film 
velocity (Equations (11-25), (11-28) and (11-48)) results in 

h. — 1 = d \h h Ar u - f(i,\ 5 l (ii-49) 

K 6 *1 1, *» m if ■ 

Using the relationship between dx and d<f> from Equations (11-14) and (11-15), and 

carrying out the differentiation results in 

A, A T i P/ A, M U m i 



h f 5 2t| ( h f x(4>) 



-($ fzW) 



+ / 2 ( 4 >) 



dh_ 



(11-50) 



Defining a dimensionless parameter, Z 2 , as 



S 2 



Z 2 = — 

B i 



(11-51) 



where 



B 2 = 



P, U. 



(11-52) 
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Equation (11-50) becomes 



^1^2 z = 4 (H-53) 

W) 2 D' / 2 ( 4 >) 



This problem is a first order ordinary differential equation which can be solved using 
a forward stepping numerical integration technique. The initial condition for the 
problem is based on symmetry which requires that (dZ->/d 4>)^ =0 = 0. Solving Equation 
(11-53) at 4> = 0 results in the initial condition 

7 _ 2 x(0) _ 2uk 2 



-2 4>=0 



D e f 2 (0) D e (l+k) 



(11-54) 



The non-dimensional film thickness, 6 , and film thickness, 6, are determined 
respectively by 



s' = ,/^mo 



(11-55) 



6 = . 

The heat-transfer coefficient for forced convection condensation can be considered in 
a different form. Combining Equations (11-40), (11-51), (11-52) and (11-55) results in 



Nu Re 



f*- 1 ' 2 = _L 



8 '(40 



, and 



(11-56) 



tin Re-" - 



* J ° 8’(4>) 



For the case where k = 1 (circular tube), / 2(4^) = 2sin <J>, dx = ad<\> and Equation (II- 
49) reduces to 

1 = 1 d [p i U. 

6 a d<J> I 



sin<t> 6 



(11-57) 



which is the same as the Shekriladze-Gomelauri [21] governing equation for a circular 
cylinder with no body forces. 



36 



D. MIXED CONVECTION CONDENSATION ON AN ELLIPTICAL CYLINDER 



(SHEKRILADZE-GOMELAURI [21] TYPE ANALYSIS) 

Consider a saturated vapor at temperature T sat , flowing downward over a 
horizontal elliptical cylinder with free stream vapor velocity, U w . The elliptical 
cylinder has eccentricity, k, and semi-major axis, a, which is parallel with the direction 
of gravity and vapor flow. The forces acting on a film element are viscous and body 
forces. The interfacial shear is given by the Shekriladze- Gomelauri [21] model for 
infinite condensation rate as in the case of the forced convection model (Equation (II- 
42)). 



Conservation of mass and energy are as previously derived. Conservation of 
momentum for the condensate film reduces to 



A, ^-7 + Pi 8 /i(4>) = 0 

dy 2 



(11-58) 



with boundary conditions, 



= 0 



(fl 



m £L 



(11-59) 



'y=b 

Integrating the momentum equation, subject to the boundary conditions, results 
in the local film velocity 



«O0 = — /,(4>) I i >y-~ 

\ 2 ) 



m U 

+ - / 2 (<i>) y 



(11-60) 



’ll 



and mean film velocity 



Pi 8 
’ll 




U m 

m 



(11-61) 



where the expression for potential flow over an ellipse and its associated velocity 
function (Equations (11-22) and (11-47)) have been used. 
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Combining the continuity equation, energy equation and expression for mean film 



velocity 



(Equations (11-25), (11-28) and (11-61)) results in 

x i 1 _ d (pfe 8 3 P, AT 
h fg 8 dx\ n, 3 t,, h fg 




(11-62) 



Using the relationship between dx and d<j>, and algebraically manipulating Equation (II- 
62) results in 



6 * 



d_\ 

2 x(<t>) 



2F/ 1 (0) 




+ m) a* 



(11-63) 



where F is a dimensionless parameter relating free and forced convection and 8 is a 
dimensionless film thickness given respectively by 

*>« h 8 



F = 



u:k, A T 



(11-64) 



8 * = 6 



tL^,±{ge 

'i i, o. 



Differentiating Equation (11-63) results in 
1 D 



8* 2 x(<J>) 



8- J 



2F/,(<j>) 6' 2 ^ + 2F/fa>) - / 2 (4>) 

<k p 3 



<f8‘ 

<f<t> 



f 2 (4>) 8 



i- 



(11-65) 



Regrouping the terms shows that the equation is a first order ordinary differential 
equation: 

{ 2F /,(<!>) 8* s + / 2 (<t>) 6*}^ ♦ 1 2F /fa) -y + /&>)} = . (» -66) 

The initial condition for this problem is again based on symmetry, which requires that 
d 6 " / d 4>=0 at 4>=0. Inserting this condition into Equation (11-66) results in a fourth order 
expression for d* 

| F& 0) S* 4 + / 2 (0) 8* J - 2 ^ = 0 , (11-67) 

j u 
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whose rool of interest is given by 



8 *( 0 ) 



-/ 2 ( 0 ) 



N 



(f 2 (0)f * 4 <2/3 F/ftO)} {2 
2 (2/3 FrfiO)} 



1/2 



x(Q) } 



( 

3k 


(Uk)\ 


16 F a 


1+jfc 


4F\- 


\[ k ) 


3 D e 


k 

/] 



( 11 - 68 ) 



Equation (11-66) can be solved numerically to obtain the film thickness. The other heat 
transfer parameters are determined as before (Equations (11-39) and (11-56)). 

For the case where F=Q and k = 1 (forced convection over a circular tube), Equation 
(11-63) reduces to Equation (11-57) as expected. For the case where F - »and k =1 (free 
convection over a circular tube), Equation (11-63) becomes 



a* 



d < d_ 

2x(4>) <*t> 



{I ^/,(4>) 6 




(11-69) 



Upon substitution for the parameters F and 6 (Equation (11-64)), Equation (11-69) 
becomes 



J. 

a 



1 P< 2 \ 8 d_ 
a A, A T d§ 




(11-70) 



which is the same as Nusselt's governing equation [5] for laminar film condensation on 
a horizontal circular cylinder. 



E. MIXED CONVECTION CONDENSATION ON AN ELLIPTICAL CYLINDER 
WITH SURFACE TENSION AND PRESSURE GRADIENT EFFECTS 
Consider a pure saturated vapor at temperature T sat , flowing downward over a 
horizontal elliptical cylinder with free stream vapor velocity, U ro . The elliptical 
cylinder has eccentricity, k , and semi-major axis, a , parallel with the direction of 
gravity and vapor flow. In addition to the assumptions used for mixed convection in 
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the last section, pressure effects in the condensate film are considered in the 
development of the governing equations. The pressure on a condensate film element 
will be a combination of that impressed upon the element by the potential flow of the 
vapor as well as that due to surface tension. A fluid element is shown in detail in 
Figure 11*5. 




Conservation of mass and energy are as previously described. Conservation of 
momentum for the film element is a balance of viscous, pressure and body forces which 
reduces to 



with boundary conditions 



d 2 u 



dy 



+ P i 8 /i(4>) 




u 



y - 0 



= o 



du] m m b 

n, 



(11-71) 



(11-72) 
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As previously stated, the film pressure is a combination of the pressure due to 
potential flow, p^, and the pressure due to surface tension, p a . Thus the total pressure 
gradient is given by 

dp = + (11-73) 

dx dx dx 



The pressure gradient due to potential flow is given by 

_ _ u dU ± _ -Pv^ ^ 
dx " Pv * dx x(4>) 

Substituting the expressions for potential flow and x(<J>) (Equation 
into Equation (II-74) results in 



(11-74) 

(11-22) and (II- 15)), 




dx 



p v u: 

2a 



m) 



(11-75) 



where 



/ 3 (<t>) = 



(1 +k) 2 k 1 sin2<t> 

[ sin 2 4> + <r 2 cos 2 4> ] 2 



cos 2 4> + fc 2 sin 2 4> 




+ —(A: 2 -!) 2 sin^ 



The pressure due to surface tension is given by 

P a = -2- . 

9 m) 

It is assumed that 6«R. Thus, the pressure gradient is given by 

dP„ a dR _ a dR 
dx R 2 dx R 2 x (<j>) d<\> 



(11-76) 



(11-77) 



(11-78) 



Substituting Equations (II - 9) and (11-15) into Equation (II -78) results in 



4P<, 3 a .... 

— = / 4 ( 4 >) 



dx 



2 a 2 



(11-79) 
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where 



/ 4 (4>) = 



k (1-A: 2 ) sin24> 

[ sin 2 4> + fc 2 cos 2 4>] 5/2 



\ 



cos 2 <{) + fc 2 sin 2 4> 
k 2 + (Jt 2 — l) 2 skftty 



(11-80) 



Integrating the momentum equation, subject to the boundary conditions results in 
a local film velocity of 



u(y) = 



P 1 8 .... P v , ... 3 a .... 

/j(4>) + / 3 (<t>) + - — - / 4 (<t>) 






2 a ti, 



2 T | .<2 2 



by 



-z! 



+ — / 2 (4>)y (II ‘ 81) 

"Hi 



and a mean film velocity, 



u = 



P/ & 3 a . ... 

/i(40 + — / 3 (<t0 + — - / 4 (4>) 






2n n, 



2 T) i a 2 



§2 m f/ 



(11-82) 



Combining the continuity equation, energy equation and expression for mean film 
velocity (Equations (11-25), (11-28) and (11-82)), results in 



Ar J. = d_ 
11 ^ 6 dx 



2 

Pi 8 

n, 



/,(<*>) + p \[ vU ~ m ) + 4 73/4(40 



2n ti, 



3 P; 0 
il, a 2 



& 

3 



(11-83) 



p, A, Ar y. 

2 hft n, 



/ 2 (4>) 8 



Using the relationship between dx and d<f>, and algebraically manipulating Equation (II- 
83) gives 



J_ = D < d 
b' X(4>) d<f> 



D ' 3 ( D 

f /j(4>) + pm ) + 4 — 

2a 2 \ a ) 



\2 



Bo 



m) 






(11-84) 



+ 4 A(40 6 ‘ 



where F and 6 are as previously defined in Equation (11-64). P and Bo are 
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dimensionless parameters defined by: 



p _ Pv h Ji 
P, AT 






Pi * ^ 



(11-85) 



P is related to the pressure gradient and Bo, the Bond number, relates the inertia effects 
to surface tension effects. Differentiating Equation (11-84) results in the following first 
order ordinary differential equation: 



2 m 

D. 



2 p /,(4>) ♦ — pm) * 3 

a 



'D\ 2 F 



- i-m) 

a I Bo 



6" 






2 f m) * p^(4>) 



( 11 - 86 ) 



, f D e ) 
+ 3 | — 
a 



Bo 



■y + ^(<1>) 



In the previous development, film separation and instability were not applicable 
since no effects were considered which could retard the film flow and (d6 /d<{>)-<» only 
at 4>c=rc. However, in the present analysis, the decreasing potential flow and increasing 
radius of curvature of the elliptical cylinder causes the pressure gradient to have a 
retarding effect on the condensate film over the back half of the tube. At the point of 
condensate film separation from the tube wall, the film shear stress at the wall is 
negative which implies flow reversal (du/dy ^ 0 at y = 0). Differentiating Equation (II- 
81) and solving for the film separation criteria (du/dy = 0 at y = 0) results in the 
following condition for separation to occur: 

2 Ff x ($) + ^ P / 3 (<t>) + 3 
a 






6* J + 2/ 2 (<|>) 8* £ 0 



(11-87) 
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Condensate film instability occurs when the rate of change of film thickness is infinite. 
From Equation (11-86), this condition occurs when 




( 11 - 88 ) 



From a comparison of Equations (11-87) and (11-88), it is evident that the condition of 
film instability will occur first (has a smaller magnitude than the film separation 
criteria for all conditions) and is therefore the limiting condition. Rose [28] suggests 
that this instability may manifest itself in wavy or turbulent flow, or may signify the 
detachment of the liquid film from the tube surface. For the case where 1/Bo = 0 
(surface tension neglected), and using Equations (11-21), (11-47) and (11-76) for / 1 (<J)),/ 2 ((t>) 
and / 3 (<J>) in Equation (11-88), it can be seen that the condition for which solutions can 
be determined over the entire tube surface occurs when 



For a circular tube (Ar= 1), Equation (11-89) reduces to P<F/ 8, which was that found by 
Rose [28]. If the opposite of Equation (11-89) is true, then the pressure gradient effect 



result in a rapid thickening of the condensate film. Note that the condition for 
condensate film instability is a function of vapor velocity ( F ) since the pressure 
gradient due to potential flow is velocity dependent. 

For the case where P = 0, and using Equation (11-88), the condition which allows 
solution of the momentum equation over the entire tube is now given by 




(11*89) 



is dominant enough to significantly retard the flow over the back half of the tube and 



1_ < \ 2 k A 




(11-90) 



Bo 3 !Z)I 1 -k 2 
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Note here that the condition for film instability is only a function of geometry. 
Therefore, for small values of eccentricity (approaching a vertical flat plate), the 
retarding effect of surface tension on the condensate film over the back half of the 
tube is significant and can result in a rapid thickening of the condensate film at 4>c* 
The initial condition for this problem is again based on symmetry, which requires 
that (dd*/d<J))^ 0 = 0. Combining this condition with Equation (11-86) results in a fourth 
order expression for 6 : 



Equation (11-86) can be solved numerically to obtain the condensate film thickness. The 
other heat transfer parameters are determined as before (Equations (11-39) and (11-56)). 

It can be shown that this formulation reduces to the models developed for circular 
tubes. Consider the case where k - 1 and 1/Bo=0 (circular tube with no surface tension); 
Equation (11-84) reduces to: 




whose root of interest is determined by 




B =/ 2 (0) 

C = 2 x(°) 



(11-92) 



W-4 AC-8 

\ 2 A 




(II -93) 
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which is identical to that developed by Rose [28] where 



6 * = 6 




v/2 



£ * 

° Rose 



(11-94) 



For the case where k= 1, P=0 and l/Bo=0 (circular tube, no surface tension and no 
pressure gradient), Equation (11-84) reduces to the Shekriladze- Gomelauri [21] model. 



F. MIXED CONVECTION CONDENSATION ON AN ELLIPTICAL CYLINDER 
WITH VAPOR BOUNDARY LAYER SEPARATION (FUJII ET AL. [22] TYPE 
ANALYSIS) 

The mixed convection model developed by Shekriladze and Gomelauri [21] has 
been used because of its simplicity and ease of solution. However, since the interfacial 
shear is approximated by an asymptotic expression based on an infinite condensation 
rate and uses potential flow outside the vapor boundary-layer (which is always 
positive), vapor boundary- layer separation does not occur. As described in Chapter I, 
Fujii et al. [22] modified the interfacial shear stress expression by simultaneously 
solving the boundary-layer equations for the condensate and vapor ensuring compatibil- 
ity at the condensate/ vapor interface. This technique allows for a more precise 
description of the condensation problem under forced convection conditions and 
includes vapor boundary-layer separation. 

Consider a pure saturated vapor at temperature T sat , flowing downward over a 
horizontal elliptical cylinder with free stream velocity, U m . The elliptical cylinder is 
oriented such that its major axis is aligned with the direction of gravity and vapor flow. 
Due to the complexities of solving the two-phase boundary-layer equation, it is assumed 
that the only forces acting on the condensate film element are gravity and viscous 
forces. Surface tension and pressure gradient are neglected which is reasonable since 
they were found to have negligible influence under most conditions (discussed in 
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Chapter IV). Other assumptions are the same as in Section II. E. The vapor is assumed 
to be in thermal equilibrium such that the energy equation for the vapor need not be 
considered in the analysis. The only forces acting on an element in the vapor boundary- 
layer (see Figure II - 6) are inertia, viscosity and pressure gradient, which is impressed 



upon the element by potential flow outside the vapor boundary- layer. The tangential 




Figure 11-6. Condensate Film and Vapor Elements for Mixed Convection using 
Fujii et al. [22] Type Model. 



velocity at the vapor/condensate interface is much smaller than the potential velocity 
outside the vapor boundary-layer and can therefore be assumed to be negligible. Based 
on the above assumptions, conservation of mass and momentum for the vapor boundary- 
layer are given by 



au 

dx 



+ *? = o 






dx 



dy 



dy 

dx 



0 continuity ) (a) 



n v #u , „ . ... 

+ — ( momentum ) ( b ) 

P„ dy 2 



(II -95) 
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For the condensate film, conservation of mass, momentum and energy are given by 



Hi 

Pi dy 2 



— — - + g /j(<|)) = 0 ( momentum ) (a) 



^1 ^sat ^wa ll> 



(11-96) 



m = ~dx ^ 1 Um = ‘ ' h* i — ( ener sy and maS5 ) (*) 

fs 



where conservation of mass and energy have been combined as in previous sections. 
Boundary conditions for this system of equations are given by 



“>=0 = VO = 0 



(fL'° 



(11-97) 



where A is the vapor boundary-layer thickness. The compatibility conditions at the 
interface require that 

u ’-‘ ° ■ -(fL.-qfL.-- 

h ^sat~^wall) 



(II -98) 



and - p v V y . 4 = 



K 6 



where the last compatibility relation is the condensation mass flux rate (i.e. for a given 
surface area, the amount of vapor condensed must equal the increase in the mass of the 
film). 

The technique for solving the vapor boundary-layer momentum equation involves 
an approximate integral technique developed by Truckenbrodt [23]. The initial problem 
development follows the method of Pohlhausen [36]. The continuity and momentum 
equations (Equations (II-95a) and (II - 95b)), are integrated over the thickness of the 
boundary layer from y = 6 to y = 6 + A and then combined to eliminate the normal 
component of vapor velocity, K, resulting in a momentum integral equation 

1 d (ww2 



j_ fL ((,’ 4,1 * £5i ■ 

ui * ' ", * u . p. u, 



(II -99) 



4> P v u * 
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where Ajand Ao are the vapor boundary layer displacement and momentum thicknesses 
given respectively by 



and 




(II- 100) 



V 5 is the normal component of vapor velocity at the vapor/condensate interface due to 
the condensation process itself and is sometimes called the vapor boundary-layer suction 
velocity. Differentiating the A-> term in Equation (11-99) and rearranging the equation 
results in: 



d*2 + A* 



dx 



U A dx 



2 + — 1 
A 2 ) 



U* 



Pv ul 



(II- 101) 



With some further algebraic manipulation this yields 



-\ 


W* Pvl 


i = — 


T J _ 


A,1 
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A 2 P v dU, V s A 2 p v ‘ 


dx] 


1 ->■ ] 


1 U * 


-=1 
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A , 


flv dx tlv 



(11-102) 



To solve this expression, the vapor velocity distribution across the boundary- layer 
thickness must be determined. This distribution is a function of a pressure gradient 
parameter, k, and a suction parameter, k 1? given by Tuckenbrodt [23] as 

A 2 Pv (11-103) 



^ P v 

K = 

Tl v dx 



and Kj = 



Truckenbrodt [23] defined a shear function f(K,Kj) and a shape function H(k,Kj) as: 

A, 



/(*, K,) = 



A 2 
Ov U, 



and /f(K,Kj) = 



(11-104) 



and a dimensionless function, Z, given by 



Z = 







(11-105) 
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Equation (11-102) becomes 



D — = 2— f / - (2+H) k - K.1 . (11-106) 

* dx U^'- ‘J 

Based on existing theoretical data (Schlicting [37] and Torda [38]), Truckenbrodt [23] 
proposed a simple linear approximation for Equation (II- 106) given by 

D — = — [0.441(1 -2k.) - 6k I . (11-107) 

dx U , L J 

Substituting Equations (II- 103), (II- 105) and (II- 106) into Equation (II- 107) and defining 
dimensionless velocities by 

U \ = ^ and V. = — , (II- 108) 

* u 5 u 



results in 



— = XGM | 0 441 f ! +2V i jRe^Z ] - 6 —L.^± Z 

d* U^D, \ L »V v j m ^ 



(11-109) 



The dimensionless suction velocity, V 6 , is determined by transforming the compatibility 
relation involving the condensate mass flux rate in Equation (11-98) and is given by 

O 



where Gis a dimensionless parameter which is proportional to the condensation rate and 
is defined by 



G = 



(^sat 



n, h A p v Tj 



p« n, 



(ii-ni) 



Combining Equations (11-109) and (II- 110), the governing equation for the vapor- 
boundary layer becomes 



dZ _ l x(4» 

^ D e 



1 0.441 


1 - 2—yfZ 
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6' 



- 6 



D, dU ' 
x( 4 >) <*<t> 




( 11 - 112 ) 
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The initial condition for Equation (11-112) comes from symmetry and is given by 

(—) =0 . (11-113) 

Un = 0 

Combining Equations (11-112) and (11-113) results in 



6 A K' 



, ns . m f + 0 882 G — (°) \^(0) - 0.441 = 0 

X(0) l 5* 



(11-114) 



This is a quadratic equation from which Z( 0) is determined to be, 



Z( 0) = 



-0.0735 G x(0) 


K) 




l d<\> j 
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/ 



0.0735 



V 



(M*) 


" l ? 


U<t> j 


*=0 8 *(0) J 






(11-115) 



The interfacial shear is determined from Equation (11-104). Truckenbrodt [23] 
provides a simple approximation for the shear function, /(k,!^), given by 



where 



/ (k,Kj) = 3.22 (K a +K) 



K a = 0 0682 + 0.174 k, 



k and Kj are modified by combining Equations (11-103) and (11-105): 



D dU. G _ 

k = — * Z W k, = — Vz . 



x(4>) # 



(11-116) 



(11-117) 



(11-118) 



5* 



Dimensionless interfacial shear is defined by: 

; _ UK- 



(11-119) 



p v u: 



An expression for dimensionless shear can then be obtained by combining Equations (II- 
104), (11-105) and (11-116): 



t/ A 



-c. = 6.44 — - Jk (k+k) 

4 fz v a 



(11-120) 
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The governing equations for the condensate film are solved in the usual manner. 



The momentum equation (Equation (II-96a)), is integrated to solve for condensate film 
velocity, 

U = /,(<D) L - l!) + h y . (II -121) 

n, v ti, 



The mean velocity is determined from Equation (11-24): 

/,(4>) 



— 52 + A 6 



h, 



2n, 



(11-122) 



Substituting Equation (11-122) into the energy/continuity equation (Equation (II - 96b)) 
and multiplying both sides of the equation by 

11 ' D< h * (11-123) 

Pf Ul h 

results in 
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(11-124) 



Carrying out the differentiation of Equation (11-124) and rearranging results in the 



ordinary differential equation, 



d y 

dt |) 



X(<t>) _ f 6* 4 - 1 1 5* 3 

D e 3 d$ 4 G d<$ 

Fm 6* 3 ♦ j 1 x s 6* 5 



(11-125) 



which, due to the symmetry of the problem, has the initial condition, 
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(11-126) 
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(11-127) 



Combining Equations (II 
X(0) 



-125) and (11-126) results in 
3 U<t>U 4 G[d* 



6 * 3 ( 0 ) = 0 



4>=0 



from which 6*(0) can be determined by numerical methods. 

For the case of k - 1 (circular tube), Equation (11-125) reduces to the analysis of 
Fujii et al. [22] as corrected by Lee and Rose [25]. The differential equation for the 
vapor boundary-layer (Equation (11-112)) and for the condensate film (Equation (II- 
125)) are solved simultaneously using the compatibility relations of Equations (11-117), 
(11-118) and (11-120). At the point of vapor boundary-layer separation, the interfacial 
shear stress becomes negative. Downstream of this point, it is assumed that the shear 
stress at the vapor/condensate interface is negligible and a simple Nusselt type analysis 
is used for the remainder of the elliptical surface. The heat transfer parameters are 
determined as before (Equations (11-39) and (11-56)). 
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III. NUMERICAL METHOD 



In the theoretical development, the governing equations for each model were 

either reduced to a single first order ordinary differential equation (ODE) or to a 

simultaneous solution of a system of first order ODEs of the form 

$=>(x,y) • (Ill - 1) 

dx 

Initial conditions were determined from the symmetry of the problem which required 
that both the slopes of the condensate film and vapor boundary- layer were zero at <J> = 
0. These systems of equations are solved using one of the forward stepping numerical 
methods. 

A. SOLUTION TECHNIQUE FOR THE NUSSELT [5] AND SHEKRILADZE- 
GOMELAURI [21] ANALYSIS METHODS 

Each analysis method resulted in a single first order ODE of the form given by 
Equation (III - 1). The solution of this equation provides the condensate film thickness 
along the elliptical tube surface. Many numerical methods are available to solve this 
type of problem, each with its own advantages and disadvantages. For this study, a 
fifth order Adams method predictor-corrector algorithm was developed from Crandall 
[39). The Adams methods are multi-step techniques which use finite difference type 
operators incorporating previously determined points (hence the term multi-step) to 
determine the value corresponding to the next step. These methods provide the same 
accuracy and are more efficient than the Runge-Kutta (single step) methods but are 
restricted to a fixed step size. The Adams method employed in this study uses an 
Adams-Bashforth explicit method [39] to predict the value of the dependent variable at 
the next step and then an Adams-Moulton implicit method [39] to correct this predicted 
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value. This method is not self starting and requires determination of the first four 
points. To start-up the algorithm, lower order Adams methods were used to obtain these 
four points. The recurrence formulae for this algorithm are given by 

y. = y*-i + h E PjMw) (in - 2 ) 

k 

where x is the independent variable, y is the dependent variable, h is the step size, s is 
the step number for which y is to be determined, k are the points determined in 
previous steps and P k are the coefficients provided in Table III -1 and III - 2. 



TABLE III- 1. RECURRENCE COEFFICIENTS FOR PREDICTOR FORMULA 



FORMULA 


Pk 

s-5 s-4 s-3 s-2 s- 1 


TRUNCA- 

TION 

ERROR 


(1) 


i 


0(h 2 ) 


(2) 


-1/2 3/2 


0(h 3 ) 


(3) 


5/12 -16/12 23/12 


o(h 4 ) 


(4) 


-9/24 37/24 -59/24 55/24 


o(h 5 ) 



TABLE III - 2. 


RECURRENCE COEFFICIENTS 


FOR CORRECTOR 


FORMULA 


FORMULA 


p k 

s-4 s-3 s-2 


s - 1 


s 


TRUNCA- 

TION 

ERROR 


(1) 




1/2 


1/2 


o(h 3 ) 


(2) 


-1/12 


8/12 


5/12 


0(h 4 ) 


(3) 


1/24 -5/24 


19/24 


9/24 


o(h 5 ) 


(4) 


-19/720 106/720 -264/720 


646/720 


251/720 


o(h 6 ) 



The algorithm used in the solution for mixed convection condensation with surface 
tension and pressure gradient (Section II.E) is provided in Appendix (B). 
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For many applications, this type of start-up routine would normally not provide 
sufficient accuracy. However, no difficulties were encountered using this numerical 
procedure for solving the asymptotic interfacial shear stress models. Solution 
convergence was obtained for an angle step size of 0.1°. The difference between an 
angle step size of 1.0° and 0.1° was less than 1%. 

B. SOLUTION TECHNIQUE FOR THE FUJII ET AL. [22] ANALYSIS METHOD 

Reduction of the governing equations resulted in a system of two simultaneous 
ODEs, one for the momentum thickness of the vapor boundary-layer and one for the 
condensate film thickness. The Adams method previously described did not provide 
sufficient accuracy to solve this two-phase boundary layer problem. In particular, the 
system of equations exhibited a stiffness problem, i.e., the numerical method could not 
accurately solve the vapor boundary-layer ODE over the first several steps even with 
an angle step size of 0.001°. To obtain sufficient accuracy, a much smaller step size 
would be required which would significantly reduce the efficiency of the computer 
algorithm. Lee [40] indicated that similar difficulties had been encountered when Lee 
and Rose [25] solved the same equations for a horizontal circular cylinder. Lee and 
Rose [25] used a Runge-Kutta method which permitted variable step size. Thus, the step 
size could be reduced sufficiently to obtain the required accuracy for the vapor 
boundary-layer ODE over the front of the tube and then increased over the remainder 
of the tube to improve the efficiency of the algorithm. Utilizing this numerical, 
procedure, they were able to solve the system of equations over most of the range of 
dimensionless parameters (F and G). However, in some cases (specifically large values 
of F and G), the system of equations were still too stiff for solution. 

To overcome the problem of stiffness, a numerical problem solver (IVPAG) from 
the International Mathematical and Statistical Library (IMSL) [41] was used. This 
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algorithm uses Gear's stiff method [42] to efficiently solve systems in which stiffness 
is a problem. However, as in the study of Lee and Rose [25], this method could still not 
solve the two-phase boundary-layer equations over the entire range of dimensionless 
parameters considered, but did cover the practical ranges of the parameters. 

It was originally conceived that the computer algorithm would progress around 
the tube surface until the criteria for vapor boundary -layer separation was reached 
(dx^/d^ < 0). At this point, the model would shift from solution of the two-phase 
boundary-layer equations to the simple Nusselt model. However, the numerical solver 
encountered difficulties in solving the system of equations as it approached the point 
of vapor boundary-layer separation. Therefore, the solution method was changed to a 
two step procedure. In the first step, the system of equations were solved until the 
computer algorithm indicated a problem had been reached (usually the algorithm 
indicated a stiffness problem had been encountered). The interfacial shear stress data 
was analyzed to verify that this problem was the result of rapidly changing conditions 
associated with vapor boundary-layer separation. Therefore, to facilitate the solution, 
vapor boundary-layer separation was assumed to have occurred if the interfacial shear 
stress at the last point obtainable was less than ten percent of the maximum interfacial 
shear stress observed over the surface. At this point, the problem solver was within one 
to two step sizes of the actual separation point (for a step size of 0.1°, the computed 
separation point was therefore within 0.2° of the actual separation point). The second 
step involved re-running the analysis using the computed separation point calculated 
in the first step as the point at which the two- phase model was switched to the simple 
Nusselt model. A converged solution was obtained for an angle step size of 0.1°. The 
algorithm for solving the two-phase model is provided in Appendix (C). 
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IV. RESULTS AND DISCUSSION 

When comparing solutions for an elliptical tube with those of a circular tube, 
equivalent surface areas have been used. The effective diameter of the equivalent 
circular cylinder is given by Equation (11-16). The dimensionless streamwise length x * , 
defined as the ratio of the streamwise length, x, to the half perimeter length and given 
by 

x* = , (IV-1) 

« D , 

is used when comparing local values of film thickness and heat- transfer coefficients. 
Its use enables direct comparison between elliptical and circular tubes since it represents 
an equivalent fraction of the total perimeter from the top of the tube to the 
dimensionless point. 

A practical range of eccentricities for an elliptical tube with major axis aligned 
vertically is 0.3 < k < 0.6 . The lower limit of 0.3 is based on discussions with a tube 
manufacturer (Reference [43]) in which the elliptical tube is formed by pressing a 
circular tube with a roller assembly to obtain the proper major and minor axis 
dimensions. Roller contact points would be evenly spaced to maintain the elliptical 
curvature and to prevent an overly flat surface. The upper limit is chosen as the 
eccentricity in which noticeable effects were observed in the heat transfer 
characteristics (as discussed in this chapter). 

A. EFFECT OF GRAVITY 

As previously noted (Section I.A), some improvement in condensation heat transfer 
is expected for an elliptical tube in quiescent vapor compared to a circular one since 
more of the surface is aligned with gravity. By increasing the "effect" of gravity on the 
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condensate film, the mean condensate velocity at a given streamwise location is 
increased which results in a thinning of the condensate film and increased heat 
transfer. The gravity function, f^cf)), defined in Equation (II -21) is shown in Figure IV- 
1 . 




As the elliptical tube approaches a flat plate configuration (k - 0), the effect of gravity 
increases to sweep the condensate film along the tube surface. The effect of gravity on 
the local film thickness is shown in Figure IV-2for two eccentricities. As can be seen, 
the condensate film is thinned over the front and rear portion of the tube compared to 
a circular tube but is slightly thicker in the middle region. The thickness of the film 
is controlled by the rate of condensation (which is dependent on condensate film 
thickness) and the condensate velocity (which is also dependent on condensate film 
thickness and gravity). Near the top of the elliptical tube, the larger gravity component 
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increases the condensate velocity relative to that of a circular tube resulting in a 
thinner condensate film. This thinner film, however, results in an increased 
condensation rate which tends to thicken the film further downstream. Since the 
elliptical tube has a greater condensation rate near the top of the tube, the film thickens 




Figure IV-2. Local Film Thickness for Varying Eccentricity ( k ). 



more rapidly compared to that on a circular tube causing the relatively thicker film in 
the middle region as seen in the figure. At 90°, the gravity component is the same as on 
a circular tube and so has no relative thinning effect. Over the rear half of the tube, 
the reduced condensation rate (as a result of this thicker condensate film) and increased 
velocity due to the larger gravity component relative to a circular tube results a 
continued thickening of the film on the elliptical tube, but at a slower rate than that 
on a circular tube. Eventually, the elliptical tube film thickness is again thinner than 
on a circular tube for the same dimensionless streamwise length. The condensate film 
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thickness is therefore determined by a balance of the effects of two effects: increased 

condensate velocity (which tends to thin the film) and increased condensation rate 

(which tends to thicken the film). 

The overall effect on the mean heat- transfer coefficient (given as a mean Nusselt 
number) is shown in Figure IV-3. To compare the elliptical tube with a circular tube, 

the leading coefficient of the Nusselt solution (Equation I - 3b) is plotted against 




eccentricity. Solutions for a circular cylinder ( k - 1.0) and a vertical flat plate ( k -* 0, 
L = 2a) agree well with the results of Nusselt [5]. For the horizontal flat plate ( k - «>), 
the heat- transfer coefficient approaches zero as would be predicted by the Nusselt 
model. For the practical range of eccentricities (0.3 < k < 0.6), it can be seen that the 
effect of placing more of the tube surface in the direction of gravity is to increase the 
mean heat- transfer coefficient by 7% for k = 0.6 and 13% for k = 0.3. Conversely, by 
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placing more of the tube surface perpendicular to the direction of gravity (i.e., 
approaching a horizontal flat plate), the mean heat- transfer coefficient decreases as 
expected. These results are in agreement with earlier elliptical tube studies [32,33]. 

B. EFFECT OF VAPOR VELOCITY 

1. Asymptotic Interfacial Shear Stress Approximation 

The streamlined shape of the elliptical tube has an effect on the vapor flow 
over the tube. The vapor velocity function, f 9 (<j>), defined by Equation (11-47) is shown 
in Figure IV-4 versus dimensionless streamwise length . As seen in the figure, the 
elliptical tube experiences higher vapor velocities than a circular tube at the front and 
rear of the tube but a lower vapor velocity in the central region. Using a balance of the 




Figure IV-4. Vapor Velocity Function, f 2(4>)** for Varying k. 
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factors which thicken and thin the condensate film, for the case of pure forced 
convection (F - 0, no gravity effect), the higher relative vapor velocity at the top of the 
elliptical tube results in a larger interfacial shear stress and hence a thinner condensate 
film and larger condensation rate compared to a circular tube. With more condensate 
flowing into the middle section of the elliptical tube and a lower relative vapor 
velocity, the condensate film thickens more rapidly than in the case of a circular tube. 
As the condensate flows over the rear of the tube, the decreased condensation rate (due 
to the thicker film in the middle region) and increased vapor shear results in thinner 
film compared to the circular tube. Figure IV-5 shows these effects in detail for a 




Figure IV-5. Effect of Vapor Velocity on Local Film Thickness for a Circular 
Tube and Elliptical Tube ( k =0.6). 



circular tube ( k = 1) and an elliptical tube ( k = 0.6) using the asymptotic interfacial 
shear stress (Shekriladze-Gomelauri) approximation and the two- phase boundary layer 
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shear stress (Fujii el al.) approximation (to be discussed later). Whereas the gravity 
effect (f ! ( 4 >)) described in Section IV. A was always larger for the elliptical tube (see 
Figure IV- 1), the vapor shear effect f 0 ( <J>) is larger or smaller depending on the section 
of the tube to be analvzed (see Figure IV-4). The overall effect of vapor shear using he 
Shekriladze- Gomelauri [21] model is to reduce (by about 2% in the practical range of 
eccentricities) the elliptical tube mean heat- transfer coefficient relative to an 
equivalent circular tube. 

In the mixed convection region, the mean heat- transfer coefficient of an elliptical 
tube is increased or decreased compared to a circular tube depending on the relative 
magnitudes of the vapor shear and gravity effects. This relative magnitude is measured 
by dimensionless parameter, F. For large F, gravity is dominant and the heat- transfer 
coefficient trend is described in Section IV. A . For small F , vapor shear is dominant 
and the trends described in this section determine the heat transfer. Figure IV-6 shows 




Figure IV-6. Effects of Vapor Velocity on Mean Nu Using the Asymptotic 
Interfacial Shear Stress Approximation Analysis Method. 
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the effect of F on the mean heat- transfer coefficient for varying eccentricity. As F -* 
corresponding to free convection condensation, solutions for elliptical and circular 
tubes match those reported in Section IV.A. For the case of a circular tube for all 
values of F, the present numerical solutions are within 0.4% of an equation by Rose [28]: 



Nu Re' m 



0.9 + 0.728 F 1/2 
(1+3.44 F 112 + F ) 1/4 



(IV-2) 



2. Two-Phase Boundary -Layer Shear Stress Approximation 

Vapor boundary-layer separation is not predicted by the asymptotic shear 
stress approximation since the interfacial shear is based on potential flow outside the 
vapor boundary-layer which is always positive. Solution of the two-phase boundary- 
layer equations for the condensate and vapor, with matched shear stress at the interface, 
allows determination of the vapor boundary-layer separation point and its effect on the 
condensation heat transfer. The analysis also introduces a dimensionless parameter G 
which is proportional to the condensation rate. Figure IV-5 compares the asymptotic 
shear stress approximation with the two-phase boundary layer shear stress approxima- 
tion for a large condensation rate (large G ) for a circular and elliptical tube ( k = 0.6). 
Over the forward and middle sections of the tube, the general trends are as described 
in Section IV.B.l for the asymptotic shear stress approximation. Over the rear of the 
tube, however, as the separation point is approached, the condensate film thickens more 
rapidly due to the reduced vapor shear effect (x - 0). It can be seen that for the 
elliptical tube, the separation point occurs further downstream than for the circular 
tube and thus this rapid thickening is delayed. By delaying separation, more of the tube 
surface has a thinner film. This serves to increase the heat-transfer coefficient for an 
elliptical tube compared to a circular tube if it were the only effect considered. 
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In total, however, the mean heat- transfer coefficient is affected by a 
combination of gravity, vapor shear and boundary- layer separation. Figure IV-7shows 
the effect of F, G and k on the vapor boundary- layer separation point (x s ). As 
previously explained (Chapter III), solutions could not be obtained for all possible F and 
G combinations (particularly large F and G). In general, as F and G increase (less vapor 
shear and higher condensation rate), or k decreases (more elliptical), the vapor 
separation point shifts toward the rear of the tube. For low condensation rates (small 
G) corresponding to low vapor suction, vapor separation occurs near the positions 
obtained for single phase separation without suction. Little change occurs in the 
separation point location as vapor velocity is varied. For large condensation rates (large 
G), the vapor boundary- layer experiences high suction which shifts the separation point 
downstream (i.e., as G increases, x s * increases). As vapor velocity is decreased 
(increasing F), the separation point shifts further downstream. As k decreases (more 
elliptical), the vapor boundary- layer separation point moves further downstream when 
compared to a circular tube under all conditions of F and G. At high F and G, the vapor 
boundary- layer separation points come together at the rear of the tube for both the 
elliptical and circular tube. 

Figure IV-8 shows the influence that the above effects (vapor shear, gravity, 
condensation rate and eccentricity) have on the mean heat- transfer coefficient. The 
solutions for k = 1 agree with the results of Fujii et al. [22] as corrected by Lee and Rose 
[25]. Similar trends were obtained for the elliptical tubes. For large values of 
condensation rate (large G), the two-phase boundary- layer shear stress analysis agrees 
closely with the asymptotic shear stress analysis method as expected since the latter 
assumes an infinite condensation rate. For decreasing k, the results for an elliptical 
tube tend to show a reduction in heat transfer at low F (compared to a circular tube) 
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VAPOR SEPARATION POINT 






Figure IV-7. Vapor Separation Point Location with Varying F, G and k. 
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Figure IV-8. Effects of Vapor Velocity and Vapor Separation on Mean Nu using 
the Two- Phase Boundary- Layer Approximation Analysis Method. 



and an increase in heat transfer at high F, as previously mentioned. The cross-over 
point from reduction to improvement occurs at increasing values of F as G decreases. 
For large values of F, the results converge to the Nusselt analysis. This shift in the 
cross-over point is attributed to the change in the vapor boundary- layer separation point 
with varying condensation and vapor velocity (varying (7 and F). As seen in Figure IV- 
7, for large condensation rate ( G = 5.5) and low vapor velocity (large F ), the large vapor 
suction causes a delay in vapor boundary-layer separation such that it occurs at 
essentially the same streamwise location for all eccentricities. As velocity increases (F 
becomes smaller) the separation point shifts forward at different rates dependent on 
eccentricity which is evidenced by the diverging curves in Figure IV-7c. Since early 
vapor boundary-layer separation results in a reduction in heat transfer due to rapid 
thickening of the condensate film, the circular tube (. k = 1) experiences more of a 
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reduction in heat transfer compared to an elliptical tube. The effect of delaying 
separation is to cause a thinning in the mean film thickness for the tube. In the 
balance, where the streamwise component of gravity and the shift in vapor boundary- 
layer separation are positive effects (i.e., both contribute to a relative thinning of the 
elliptical tube condensate film compared to a circular tube) and vapor shear is a 
negative effect (i.e., lower shear contributes to a relative thickening of the elliptical 
tube condensate film compared to a circular tube), delays in vapor separation move the 
cross-over point from gravity dominant to shear dominant heat transfer. For the case 
of low Gin Figure IV-7a, there is no noticeable shifting of the vapor boundary-layer 
separation point as F is varied. Thus the vapor boundary-layer separation point 
contribution does not change with vapor velocity ( F ) to the balance and the cross-over 
point occurs at lower velocities (larger F). The overall effect of vapor shear using this 
analysis is typically a slight reduction (< 2%) in the heat- transfer coefficient for an 
elliptical tube compared to an equivalent surface area circular tube, dependent on the 
magnitude of the condensation rate parameter, G. 

Figures IV-9and IV-lOshow the effects of vapor boundary-layer separation 
on the local film thickness and heat-transfer coefficient for a given dimensionless 
streamwise distance. For low F and G, the separation point occurs relatively early, but 
is significantly delayed when using an elliptical tube. As F increases (vapor velocity 
decreases), the separation points remain at their same respective streamwise locations. 
The areas under the curves in Figure IV- 10 represent the total heat transfer. Bearing 
this in mind, the comparison between elliptical and circular tubes is made a little easier. 
At high F (for given G), the area under the curve for an elliptical tube exceeds that for 
a circular tube while at low F , the opposite is true. However, for a given eccentricity, 
the values of x s are the same for varying F (as noted earlier). 
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As G increases (increasing condensation rate and suction), vapor boundary- layer 
separation is further delayed (compared to low G) and is virtually eliminated at high 
F. Separation is always delayed for an elliptical tube compared to a circular tube. 
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LOCAL HLK THICKNESS 
F-O.OJ, G-0.1 




LOCAL FILM THICKNESS 
F-10.0, G-0.1 




Figure IV-9a and b. Dimensionless Local Film Thickness for Varying G> F and k 
from the Two- Phase Boundary Layer Shear Stress Method of Analysis. 
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LOCRL riLn THICKNESS 
r -C.32 , G-S.3 




LOCRL HLK THICKNESS 
r-10.0, 5-5. £ 




0.0 0.3 0.4 0.6 0.8 J.C 

Din. x 



Figure IV-9cand d. Dimensionless Local Film Thickness for Varying G, F and k 
from the Two-Phase Boundary Layer Shear Stress Method of Analysis. 
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LOCRL HU/$OR7fRC) 
f-0.01, G-O.l 




LOCflL NU/SORTIRC) 
r- 10.0, C-O.J 




Figure 10 a and b. Local Nu for Varying G, F and k from the Two-Phase Boundary 
Layer Shear Stress Analysis Method. 
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10CRI NU/SORTfRC) 
r-o.01, G-5.5 




LOCRL NU/50R7fRC» 
r-IO.O, G-5.5 




Figure 10 c and d. Local Nu for Varying G, F and k from the Two-Phase Boundary 
Layer Shear Stress Analysis Method. 
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c. 



EFFECT OF SURFACE TENSION 



As noted in Section II. E, a condition exists for which d<5*/d4> becomes infinite at 
some critical angle 4> c , What is actually happening can be discerned from Figure IV- 11 
which shows the surface tension function f 4 (<j>) versus x for varying eccentricity. 
Surface tension causes a favorable pressure gradient over the front half of the elliptical 
tube and an adverse pressure gradient over the back half. The severity of the pressure 
gradient is localized to small regions at the top and bottom of the elliptical tube where 



a 




Figure IV- 11. Surface Tension Function, f 4 (<J>) for Varying Eccentricity. 



change in surface curvature is most severe. Therefore, it is most significant for small 
values of eccentricity. The above analysis explains the reason for condensate flow 
instability discussed in Section II. E and predicted by Equation (11-90). 
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Further discussion of surface tension is restricted to those cases in which a 



solution could be obtained over at least 90% of the tube surface. This restriction was 
arbitrarily chosen based on the small amount of heat transfer which occurs on the lower 
part of the tube and would therefore have minimal impact on the accuracy of the 
results. Figure IV-12 shows the effects of surface tension on the mean heat- transfer 
coefficient for k — 0.6 and 0.4 . The Bond number, Bo , gives the relative effect of 
inertia to surface tension. Values of 1 IBo of 0.01 (typical for steam) and 0.001 (typical 
for refrigerants) as well as 1/Bo = 0 have been shown on the figure. A number of 
interesting points can be highlighted. Firstly, the surface tension effect is much smaller 
for highly wetting refrigerants, as expected. Secondly, inclusion of surface tension 
leads to a small decrease in the mean heat- transfer coefficient (< 2%) over the whole 
range of F for the practical range of eccentricities, suggesting that any thinning of the 
condensate film over the top half of the tube is more than offset by a thickening over 
the lower half. As eccentricity decreases, this discrepancy is accentuated. Finally, 
surface tension effects are felt more in the free convection region (high F) than in the 
forced convection region (low F) as film thickness becomes dominated by vapor shear. 

Krupiczka [29] used the curvature of the condensate film surface to analyze 
surface tension on a horizontal circular cylinder. He determined that the effect of 
surface tension was only significant over the bottom portion of the tube. This fact 
implies that since the film is very thin over the top of the tube, its curvature is very 
close to that of the tube surface itself, which, for a circular tube results in no surface 
tension effect (constant radius of curvature). As the condensate flows around the lower 
half of the circular tube, the film curvature no longer follows the tube surface. Rather, 
the radius of curvature is increasing resulting in a small improvement in the heat- 
transfer coefficient. A similar analysis can be applied to the elliptical tube. Over the 
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Figure IV- 12. Effect of Surface Tension on the Mean Nu for k = 0.4 and 0.6 . 
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top half of the tube where the condensate film is thinnest, modeling the surface tension 
effect using the curvature of the tube surface is valid. However, the adverse pressure 
gradient over the lower half of the tube is probably not as severe as seen in Figure IV- 
11. Over this portion of the tube, the film is thicker and its decreasing surface 
curvature is not as severe as the tube surface curvature. The conclusion drawn from 
this analysis is that surface tension does not have a significant effect on the mean heat- 
transfer coefficient for an elliptical tube. 

D. EFFECT OF PRESSURE GRADIENT 

The effect of pressure gradient was analyzed using the asymptotic shear stress 
approximation (similar to Rose [28] for a circular tube). In this analysis, the pressure 
gradient in the condensate film was assumed to be due to the pressure gradient 
impressed on the condensate by vapor potential flow. As in the case of surface tension, 
there is a relationship between parameters P and F (Equation 11-89) for which a solution 
could not be obtained for the entire tube. Figure IV- 13 shows the pressure gradient 
function versus x* for varying k which sets up a favorable pressure gradient over the 
top of the tube and an adverse pressure gradient over the lower half of the tube. The 
effect of eccentricity is to shift the point of maximum pressure gradient to the front 
and rear of the tube. Thus, <J> C shifts downstream for an elliptical tube compared to a 
circular tube for similar P and F . This is a result of the more streamlined shape and the 
smaller resultant pressure drop over the streamwise length of an elliptical tube. 

As in the analysis of surface tension, discussion is limited to those cases for which 
solutions could be obtained over the entire tube. When pressure gradient was included 
in the momentum equation of the condensate film, solutions agreed with the results of 
Rose [28] for a circular tube, yielding a maximum 5% decrease in the heat- transfer 
coefficient over the entire range of F compared to the case with pressure gradient 
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Figure IV- 13. Pressure Gradient Function, f 3(4>) ^ for Varying Eccentricity. 



neglected. Figure IV-14 shows the mean heat- transfer coefficient for a circular tube 
(k = 1) for varying F and P. The favorable pressure gradient increases the mean 
condensate velocity resulting in a thinner film and an increased mean heat- transfer 
coefficient over the top of the tube. This increase is offset by the adverse pressure 
gradient which slows the condensate film velocity over the lower half of the tube, 
thickens the film and decreases the mean heat-transfer coefficient. For an eccentricity 
of 0.6, Figure IV- 15 compares the effect of pressure gradient for an elliptical tube with 
those for a circular tube for varying F and P. For both the circular and elliptical tubes, 
the effect of pressure gradient provides a slight reduction in condensation heat transfer. 
Figure IV- 15 indicates that the reduction in heat transfer was less for the elliptical tube 
compared to the circular tube as a result of the more favorable potential velocity 
distribution around the ellipse. 
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Figure IV- 14. Effect of Pressure Gradient on the Mean Nu for a Circular Tube 
(*=!)■ 



m 




Figure IV- 15. Effect of Pressure Gradient on Mean Nu, k = 0.6, Relative to a 
Circular Cylinder. 
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EFFECT OF VAPOR PRESSURE DROP 



E. 

The previous results indicate that vapor shear decreases the heat- transfer 
coefficient for an elliptical tube when compared to a circular tube. However, a 
consideration not yet discussed is that the vapor pressure drop across an elliptical tube 
is significantly smaller than a circular tube for the same free-stream velocity. 
Alternatively, for a given pressure drop, the free-stream velocity is greater for the 
elliptical tube. Since very little data is available for two-phase drag coefficients around 
tubes, single-phase data was used where the drag coefficient for the tube is due entirely 
to form drag. Values of C D for a circular (k = 1) and an elliptical tube ( k = 0.5) were 
estimated to be approximately 1.2 and 0.6 respectively (White [44]). For steam with T sat = 
60°C condensing on a horizontal tube with T wall = 40°C and U a = 25 m/s for the circular 
tube, the corresponding vapor velocity for the elliptical tube, for the same pressure 
drop, was calculated to be U M = 35 m/s. The resulting values of F are therefore 0.0257 
and 0.0132 for the circular and elliptical tubes respectively. Using the asymptotic shear 
stress assumption for both values of F gives an increase in the heat transfer for an 
elliptical tube of 16.3% when compared to a circular tube. Using the two-phase 
boundary-layer shear stress assumption gives a corresponding increase of 17.1%. 
Therefore, though the effect of vapor shear alone results in a reduction of the heat- 
transfer coefficient (as discussed in Section IV. B), when the streamlined shape of the 
elliptical tube and its effect on the vapor pressure drop is taken into account (allowing 
for higher U m for a given pressure drop), there is an increase in the heat- transfer 
coefficient when compared to a circular cylinder. 
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F. 



INSIDE HEAT TRANSFER COEFFICIENT 



So far only the outside heat-transfer coefficient has been discussed. In real 

condensers, it is the overall heat- transfer coefficient which is important and so 
consideration must also be given to the single- phase convective heat transfer occurring 
inside the elliptical tube. A preliminary survey of the literature (Incropera and Dewitt 
[45]) shows that for turbulent flow conditions where the Prandtl number is greater than 
0.5 (which is true for most practical condensers), correlations developed for circular 
tubes may be used with good accuracy to approximate the inside heat-transfer 
coefficient for an elliptical tube if the hydraulic diameter is used in place of the 
circular diameter. The hydraulic diameter is given by 



D h = 



4 A, 



(IV-3) 



P 



where A c is the cross-sectional area of the tube and pis the inside perimeter of the tube. 

The perimeter for an elliptical tube and an equivalent surface area circular tube is the 

same. The cross-sectional area of an elliptical tube is given by 

A c = n a b (IV-4) 

and is smaller than an equivalent surface area circular tube. Inside heat transfer 
correlations for turbulent flow heat transfer inside a circular tube are typically of the 
form [45] 
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(IV-5) 



Thus, the inside heat- transfer coefficient is related to the hydraulic diameter by 



a * 



- 0.2 



(IV-6) 
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As an example, consider an elliptical tube with k - 0.5 and a major axis with a = 0.01 
m. The equivalent circular tube would have a diameter of 0.0154 m. The ratio of the 
heat-transfer coefficient of the elliptical tube to the circular tube is given by 



^ elliptical _ 


elliptical 


- 0.2 


^ ^elliptical 


^ circle 


circle 




circle 



(IV-7) 



thus, 



** elliptical 
a circle 



ah 



- 0.2 



(DJ 2) 2 



= 1.035 



(IV -8) 



For this simple analysis, it can be seen that the inside heat- transfer coefficient of the 
elliptical tube is 3.5% greater than the equivalent surface area circular tube. This 
result, when added to the increase in the outside heat- transfer coefficient (for a given 
pressure drop), would increase the overall heat- transfer coefficient of the condenser. 
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V. CONCLUSIONS 



Analyses of laminar film condensation on a horizontal elliptical tube were 
conducted under conditions of free and forced convection, the latter using both an 
asymptotic interfacial shear stress approximation as well as solution of the two- phase 
boundary-layer equations. For the asymptotic shear stress approximation, the effects 
of surface tension and pressure gradient were also included. Where possible, these 
analyses have been validated against existing solutions for laminar film condensation 
on a horizontal circular tube and a vertical flat plate. 

Whether or not condensation heat transfer is improved for an elliptical tube 
compared to an equivalent surface area circular tube is dependent on a balance of the 
effects of gravity, vapor shear and vapor boundary- layer separation and the influence 
these have on the thickness of the condensate film. Under quiescent vapor conditions 
(no vapor shear or boundary- layer separation), gravity causes an increase (-10%) in the 
heat transfer on an elliptical tube compared to a circular tube of the same surface area 
due to an increase in "effective" gravity. Under conditions of forced convection, both 
shear stress approximations indicate a small decrease (< 2%) in performance due to the 
reduction in interfacial shear as a result of the better streamlined shape of the elliptical 
tube when compared to the circular tube in the same free stream vapor flow. However, 
when pressure drop effects are also considered, the higher allowable vapor velocity over 
an elliptical tube (for the same pressure drop as a circular tube) results in an increase 
in the mean heat-transfer coefficient of 15-20%. For conditions of mixed convection, 
the condensate film thickness is controlled by both gravity and vapor shear effects. The 
cross-over between heat transfer improvement (gravity dominant flow) and heat 
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transfer reduction (vapor shear dominant flow) for the elliptical tube is dependent on 
the vapor boundary- layer separation point. For large condensation rates (high vapor 
suction), vapor boundary- layer separation is delayed and the cross-over point shifts to 
higher velocities (smaller F). For low condensation rates (low vapor suction), movement 
of the vapor boundary-layer separation point is minimal, resulting in a minor effect on 
heat transfer and a cross-over point that shifts toward lower velocities (large F). 
Pressure gradient and surface tension each lead to a small decrease in the mean heat- 
transfer coefficient (< 2%) for an elliptical tube. 

In general, therefore, the outside heat-transfer coefficient of a condenser tube is 
improved by using an elliptical tube geometry. Approximation of the single-phase 
inside heat-transfer coefficient using a hydraulic diameter and a Seider-Tate type 
correlation indicates that the elliptical tube also has better inside heat transfer 
performance compared to an equivalent surface area circular tube. Thus, the overall 
heat-transfer coefficient should be improved. 
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VI. RECOMMENDATIONS 



This study represents an initial step in determining if a shell-and- tube condenser 
using horizontal elliptical tubes has better heat transfer than a comparable circular tube 
bundle. The effects of multi-tube arrangements (i.e., the influence of adjacent tubes 
and the spacing between tubes) on vapor flow characteristics needs to be evaluated to 
determine if further improvements in the overall heat- transfer coefficient can be 
achieved. It is suspected that the interfacial shear of the condensate film on elliptical 
tubes will be further reduced but the corresponding reduction in vapor pressure drop 
across the condenser will result in a net gain in heat transfer compared to a circular 
tube bundle. An analysis similar to that conducted by Aoune and Burnside [46, 47] can 
be used to evaluate multi-tube arrangements. 

Further studies should be conducted on a single elliptical tube to better determine 
the effect of surface tension on heat transfer. The surface tension model used in this 
study was based only on the curvature of the tube surface and does not take into 
account the relatively thicker film over the lower half of the tube. The curvature of 
the actual condensate film surface should be analyzed in a manner similar to Krupiczka 
[29]. This analysis, together with the asymptotic shear stress assumption, would involve 
the solution of a second order ODE. 

Analysis of the effect of offsetting the elliptical tube at an angle, a, with respect 
to vapor flow and gravity must be conducted to determine its effect on the heat- 
transfer performance of the tube. This is important since in real condensers, vapor will 
approach the tube at various angles. Additionally, since the elliptical tube is a 
streamlined body, vapor flow in a direction other than the direction of the major axis 
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may result in a lifting force on the tube which may result in tube vibration and wear. 
These effects would have to be studied and methods to minimize any resultant vibration 
must found. The problems associated with joining elliptical tubes to the condenser tube 
sheet and baffle plates also needs to be addressed. 

Finally, experiments must be conducted on elliptical tubes to validate the 
predictions of these theoretical models. 
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APPENDIX A - DERIVATION OF TANGENTIAL VAPOR VELOCITY OVER AN 



ELLIPTICAL TUBE FROM POTENTIAL FLOW THEORY 



The basic technique used to determine fluid flow velocity over a body is to 
employ a conformal transformation to translate the body from one complex plane into 
another complex plane. From this, the velocity of the transformed body can be readily 
computed using potential theory. For this particular application, the fluid velocity 
along the surface of a circular cylinder is determined and then transformed into the 
complex plane in which the cylinder has an elliptical shape. 

Consider a circular cylinder in the complex plane, z 1? as shown in Figure A-l. The 
surface of the circle is defined by 



z x = c e* e (a) 



(A-l) 



xl + y\ = C 2 (b) 

The conformal transform to translate the circular body from the z x plane to an elliptical 
body in the z 2 plane is given by 



*2 S *1 



*1 



(A -2) 



thus 



*2 - (*1 + ty) - 



1 

*1 + '>1 




+ 




(A -3) 



The real and imaginary coordinates in the z 2 plane are given by 



*2 = *i | 1 ' 



and y 2 = y t 



c 2 



(A -4) 
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Figure A-l. Circular Cylinder in z l Plane, Elliptical Cylinder in Plane- 



Substituting Equation (A-4) into Equation (A- lb) results in 






-5— -i 

HI 



(A-5) 



which is the expression for an ellipse in cartesian coordinates where the semi-major 
axis, a , is given by c + (1/c) and the semi-minor axis, b, by c - (1/c). The complex 
potential function in the z l plane for the circle is given by 



*2 J 2 



U m c* e 

h<Zj) = U m z x e * + = $ + 



(A-6) 



where 4> is the velocity field potential function and i \f is the stream function. The fluid 
velocity in the z x plane is determined from the complex derivative, 



dw „ -'T U m c 2 e 2 



f- « u . e 2 



u i ~ ,v t > 



(A -7) 
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where and v l are the fluid velocity components in cartesian coordinates. The fluid 
velocity in the z 2 plane (for the ellipse) is found by the chain rule, thus, 

dw dwdz j ( o\ 

— = = m - - iv~ . 

dzy dz x dz2 



Combining Equations (A- la),(A-2) and (A-7) results in 

^ U„ c 2 (e l <20 * n/2) - e**) 
c 2 € * 20 + 1 



(A-9) 



The magnitude of the fluid velocity is found by multiplying the complex velocity by 

its conjugate. Using Euler’s formula, 

e lQ = cos0 + i sin0 , (A- 10) 

and with some algebraic manipulation, the potential velocity of the fluid at the surface 
of the ellipse is given by 



jji _ dw dw _ c * (2 + 2cos20) 
0 dz^dz^ c 4 + 1 + 2c 2 cos20 



(A - 11) 



where the bar denotes the complex conjugate. Equation (A- 11) is further simplified 
using the trigonometric identities 

cos20 = cos 2 0 - sin 2 0 
cos20 = 2cos 2 0 - 1 



(A- 12) 



and the relationship between c, a and b, 

4 c 2 = (a+b) 2 

Thus the potential velocity is given by 

(1 +£) cos0 



U a = U 



\/cos 2 0 + fc 2 sin 2 0 



(A- 13) 



(A- 14) 



The angle, 0, is measured from the x { (horizontal) axis. The relationship between 0 and 
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4>, the angle measured from the y x (vertical) axis, is given by 

8 = — — <t> (A- 15) 

2 

Note that <{) is an angle defined in the z x plane for the circle and is the parametric angle 

defined in Section II. A . The potential velocity is 

U = v (1 + *) sine}) = v 
\/sin 2 4> + k 2 cos 2 <|> 



thus given by 

( 1 +*) 

+ k 2 cot?$ 



(A- 16) 
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APPENDIX B - COMPUTER CODE FOR ASYMPTOTIC SHEAR STRESS MODEL 



* PROGRAM MIXED1 

XXXXXXX*X*XXXXXXXXXXXX*XXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX*XX 

X 

X THIS PROGRAM NUMERICALLY SOLVES THE PROBLEM OF MIXED 

* CONVECTION CONDENSATION WITH SURFACE TENSION ON AN ELLIPTICAL 

x CYLINDER. THE PROBLEM FORMULATION USES THE ASYMPTOTIC SHEAR 
x STRESS APPROXIMATION ANALYSIS METHOD. 

X 

X THE PROGRAM USES AN AUTHOR DEVELOPED FORWARD STEPPING 
X ALGORITHM BASED ON A 4TH ORDER ADAMS PREDICTOR-CORRECTOR 

X METHOD TO SOLVE A 1ST ORDER ORDINARY DIFFERENTIAL EQUATION. 

x 

x THIS PROGRAM CAN BE RUN WITH DIMENSIONLESS PARAMETERS 
X CF,G,P, BO) SPECIFIED OR WHERE THE PARAMETERS ARE CALCULATED 

X BASED ON SPECIFIED FLUID PROPERTIES. MODIFICATIONS TO THE 
x PROGRAM (REMOVAL OF UNNECCESSARY STEPS BY "COMMENTING" THE 

x STATEMENT OUT OF THE PROGRAM) IS REQUIRED. SPECIFIED FLUID 

x PROPERTIES INCLUDE T SAT , T WALL , AND UINF. 

X 

X NUMERICAL RESULTS ARE DIRECTED TO AN EXTERNAL DATA FILE. 

X 

x MAJOR VARIABLES USED ARE: 

X 



X 


THETA 


PARAMETRIC ANGLE 


X 


X 


STREAMWISE LENGTH 


X 


DEFF 


EFFECTIVE DIAMETER 


X 


DEL 


CONDENSATE FILM THICKNESS 


X 


DELND, Z 


DIMENSIONLESS FILM THICKNESS 


X 


F,P,BO, 


DIMENSIONLESS PARAMETERS 


X 


JARPRL 




X 


UINF 


VAPOR FREE STREAM VELOCITY 


X 


TS 


VAPOR SATURATION TEMPERATURE 


X 


TW 


WALL SURFACE TEMPERATURE 


X 


DENSL 


CONDENSATE DENSITY 


X 


DENSV 


VAPOR DENSITY 


X 


CONDL 


CONDENSATE THERMAL CONDUCTIVITY 


X 


VISL 


CONDENSATE VISCOSITY 


X 


HFG 


LATENT HEAT OF VAPORIZATION 


X 


NU 


NUSSELT NUMBER 


X 


NUAVG 


MEAN NU NUMBER 


X 


ALPHA 


HEAT TRANSFER COEFFICIENT 


X 


QFLUX 


HEAT FLUX 


X 


A 


SEMI-MAJOR AXIS 


X 


K 


ECCENTRICITY 


X 


XFN 


PARAMETRIC RADIUS 


X 


G 


GRAVITY FUNCTION 


X 


FI 


VELOCITY FUNCTION 


X 


FI 1 


DERIVATIVE/ VELOCITY FUNCTION 


X 


F2 


PRESSURE GRADIENT FUNCTION 


X 


F3 


SURFACE TENSION FUNCTION 



X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

PARAMETER (PI =3. 141592654, GRAV=9. 807) 

REAL Z( 0 : 180000 ), THETA( 0 : 180000), X(0: 180000 ),H,DEFF, XI, X2,X3, 
+NUAVG,NU,NUORE,UINF,TS,TW,TL, DENSL , VISL ,TG,C0NDL,HFG,B2,DELND, 

+ DEL , QFLUX, A , K, RE, SUM, XFN , ALPHA, ALP HA 0,DQFLUX,F, PHI, 
+A1,A2,A3,A4,A5,A6,A7,A8,A9,A10,A11,TF, PSAT ,TB,TC,TD,R,TE,DENSV, 
+SIGML,FCN,P,C1,C2,C3,FC0,FC1,FC2,FC3,ZP1, ZP2, ZP3 , ZP4 , BO, 
+C4,C5,F21,F31,NUC0EFF 

INTEGER N,J,JJ,IC 

COMMON /GEOM/ A , K, DEFF/PARAM/F, BO, P 

EXTERNAL FCN, XFN , G, G1 , FI , FI 1 , F2 , F3 

0PENC20,FILE='/MIXED1 OUTPUT Al' ,STATUS='OLD' ) 
0PEN(30,FILE='/DEL2 DATA Al ' , STATUS 1 1 OLD 1 ) 
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FILE: MIXED1 



FORTRAN A 



*** INPUT FLUID PARAMETERS *** 

PRINTS, 'ELLIPSE OF THE FORM ( X/B )**2+(Y/A)**2=l ' 

PRINTS, 'INPUT "A" DIMENSION (M)T' 

READ*, A 

PRINTS, 'INPUT T SAT AND T WALL (K)?' 

READ*, TS, TW 

PRINT*, 'INPUT ANGLE STEP SIZE (DEGREES)?' 

READ*# H 

PRINT*, 'STEP SIZE, H = ' , H 
C PRINT*, 'INPUT UINF?' 

C READ*, UINF 

PRINT*, 'INPUT K? * 

READ*, K 

PRINT*, 'INPUT F?» 

READ*, F 

N= 18Q/H+1 
H=PI/180 . *H 
B0 = 0 . 

P = 0. 

*** DETERMINE FLUID PROPERTIES ******** 

Al = 15 . 49217901 
A2=-5. 6733717693 
A3=l. 4597584637 
A4=13. 877000608 
A5=-80 .887673591 
A6=123. 56333463 
A7=-188. 321212064 
A3=660. 91763485 
A9=-1382. 4740091 
A10 = 1300. 1040184 
All=-449. 39571976 
TF=T$/1000 . 

PSAT= 1 . E+6 *EXP ( Al+A2/TF+A3*L0G(TF)+TF*A4+A5*TF**2+A6*TF**3+A7* 
+TF**4+A3*TF**5+A9*TF**6+A10*TF**7+A11*TF**8) 

TB = 150 0 . /TS 

TC=2.5*LQG(1 .-EXP(-TB)) 

TD=. 00 15/(1 .+.00Q1*TS)-.000942*EXP(TB+TC)*SQRT(1./TB)-.0004382*TB 
R = 461 .51 

TE=2.*PSAT/(R*TS) 

DENSV=TE/ ( 1 . +SQRT ( 1 . +2 . *TD*TE) ) 

TL=(2.*TW+TS)/3. 

DENSL=1 ./( .0012674-TL*(2. 02915E-6-TLX3 .8333 E- 9)) 
VI$L=2.414E“5*10**(247.8/(TL-140.)) 

TG=TL/273.15 

CONDL=-.92407+TG*(2.3395-TG*(1.3007-TG*( . 52577-TG* . 07344) ) ) 
HFG=3468920 . -TS*(57 07 . 4-TS*( 11 .5562-TS*. 0133103)) 

SIGML=(-. 0003*(TL-273.15)**2-.138*(TL-273.15)+75.6)/1000. 

*** DETERMINE THETA( I ) , X( I ) , DEFF ******** 

xxx USE SIMPSON'S RULE TO NUMERICALLY INTEGRATE ALONG SURFACE *** 

THETA ( 0 ) = 0 . 

X(0)=0. 

DO 10, 1=1, N 

THETA ( I )= THETA ( 1-1 )+H 
X1=A*XFN(THETA(I-D) 

X2=A*XFN((THETA(I-1)+THETA(I))*.5) 

X3=A*XFN( THETA( I ) ) 

X(I)=X(I-l)+H/6 .*(X1+4.*X2+X3) 

10 CONTINUE 

DEFF=2.*X(N)/PI 

*** DETERMINE DIMENSIONLESS PARAMETERS AS REQUIRED **** 

C F=VISL*DEFF*HFG*GRAV/(UINF**2*CGNDL*(TS-TW)) 

UINF=SQRT(VISL*HFG*GRAV/(F*CONDL*(TS-TW))) 
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FILE: MIXED1 FORTRAN A 



C 

C 



XXX 



XXX 

XXX 



50 

60 

XXX 



P=DENSV*HFG*VISL/(DENSL*CONDL*(TS-TW)) 

30=SIGML/(DENSLXGRAVXDEFFXX2) 

RE=DENSLXDEFF*UINF/VISL 

PRINT*, 'COMPLETED THETA , X, DEFF * 

WRITEC 20 , 5000 ) 'SAT. TEMP. (K) = STS, 'WALL TEMP. = ',TW, 

+ 'I/BO = ',BO,'P = ' , P, ' F = * , F, * UINF (M/S) = ',UINF,'REL = ' , RE 
WRI TEC 20 , 5 010) * "A" LENGTH CM) = ',A, f "K f ' ECCENTRICITY = ',K, 

+ ' DEFF (M) = ' , DEFF 
WRITEC 20 , 5020 ) 

DETERMINE INITIAL CONDITION, ZO xxxxxxxx 

B2=$QRT(UINF*DENSL/(DEFF*VISD) 

THETAC 0) = 1 . E-9 

XC 0) =1 . E-9 

CALL DRV C F2, 0 . , F2I ) 

CALL DRV C F3 , 0 . , F31 ) 

CL=(2*F*GL(0. )+DEFF/AXP*F21+3*(DEFF/A)X*2*F*B0*F31)/3. 

C2 = FI ICO.) 

C3=-2.*A*XFN(0. )/DEFF 

Z(0)=SQRT((SQRT(C2XX2-4XC1XC3)-C2)/(2.XC1)) 

ALPHAO=CONDLXB2/ZCO) 

SUM=0. 

DETERMINE DIMENSIONLESS FILM THICKNESS, Z xxxx 
ADAMS-MOULTON PREDICTOR/CORRECTOR ALGORITHM XXX* 

FC0=FCN(THETA(0),Z(0)) 

ZP1 =Z( 0 )+HxFCO 

Z(1)=Z(0 )+H/2 . XC FCO+FCNC THETAC I ) , ZPI ) ) 

FC1=FCN(THETA(1),Z(1) ) 

ZP2=Z(1)+H/2.X(3.XFC1-FC0) 

ZC2)=Z(l)+H/12.XC5XFCNCTHETAC2),ZP2)+3 .XFC1-FC0) 

FC2=FCN( THETAC 2) ,Z(2) ) 

ZP3=Z(2)+H/12.x(23.XFC2-16 .XFC1+5.XFC0) 

Z(3)=Z(2)+H/24.X(9 . XFCNC THETAC 3) ,ZP3)+I9 . XFC2-5 . XFC1+FC0) 
FC3=FCN(THETA(3),Z(3)) 

DO 50, 1=4, N 

ZP4=Z(I-l)+H/24.X(55.XFC3-59.XFC2+37 . XFC1-9 . XFCO ) 
Z(I)=Z(I-l)+H/720.x(251 . XFCNCTHETAC I ) , ZP4)+646 . XFC3-264 . X 
+ FC2+106 .XFC1-19 .XFCO) 

FC0=FC1 
FC1 =FC2 
FC2-FC3 

FC3=FCN( THETAC I ) , ZC I ) ) 

C4=(2.XFXG(THETA(I) ) + DEFF/AXPXF2C THETAC I) ) + 3 . X( DEFF/A) XX2XF 
+ XB0XF3CTHETACI) ) )XZ(I)XX3+F1(THETA(I) )XZCI) 
C5=C4+F1CTHETA(I))XZCI) 

IFCC4 . LE. 0 . ) THEN 

PRINTX, 'C4 LIMIT' 

IC = I 

GO TO 60 

ENDIF 

IFCC5.LE.0. ) THEN 

PRINTX, *C5 LIMIT' 

IC = I 

GO TO 60 

ENDIF 

CONTINUE 

IC=N-1 

PRINTX, 'COMPLETED Z DETERMINATION' 

DETERMINE HEAT TRANSFER PROPERTIES XXXXXXXXX 



J = 0 
JJ = 0 

DO 100,1=1, IC 
DELND=ZC I ) 
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FILE: MIXED1 FORTRAN A 



DEL=DELND/B2 

ALPHA=CONDL/DEL 

SUM=SUM+C ALPHA OXAXXFNC THETAC I -1 ) ) +ALPHAXAXXFNC THETAC I ) ))XHX . 5 

ALPHAO=ALPHA 

IFCJJ.EQ.J) THEN 

QFLUX=ALPHAXCTS-TW) 

NU=DEFF/DEL 
DQFLUX=1 . /DELND 

WRITEC 20 , 5030) THETAC I )*180 ./PI , X( I )X2 ./ ( DEFFxp I ) , 

+ DELX1000 . , QFLUX/1000 . , DELND, DQFLUX, NU 

WRITEC 30 , X ) XCDX2./CDEFFXPI) , DELND 
J=J+1./180.XN+1 

ENDIF 

JJ=I+I 

100 CONTINUE 

WRITEC 20, 5030) THETAC 10X180 ./PI, XC IC)x2 ./C DEFFxPI) , DELX1 00 0 . , 
+QFLUX/ 1000 ., DELND, DQFLUX, DEFF/DEL 
WRITEC 3 0 , x) XCIOX2./C DEFFXPI), DELND 

NUAVG=SUMX2/CPIXC0NDL) 

WRITEC 20 , 5040 ) 'NU/SQRTCRE) = 1 , NUAVG/SQRTC RE) , * NUAVG = ' , NUAVG 

5000 FORMAT C1X,T5,70C 'x 1 ) // IX, T5, * HEAT TRANSFER PROPERTIES ' // 2C1X, 
+T1 0 , A17 , F6 . 2 / ),4C1X,T10,A17,E9.3 /)) 

5010 FORMAT C / IX, T5, » ELLIPSE GEOMETRY ' // 3 C IX, T10 , A19 , F6 . A / )) 

5020 FORMAT C / IX, T7, 'THETA 1 , T15, f DIM X 1 , T25 , f DEL f , T32 , f HT FLUX', 

+T41 , 'DIM DEL',T52, 'DIM',T58, 'NUSS #» / T7 , ' C DEG) ' , T25, ' CMM) ' , T32, 
+ * (KW/M2) ' , T50 , *HT FLUX' / T6 , 7 C ' = ' ) , T15, 6 C ' = » ) , T23 , 7 C ' = » ) , T32, 

+7C ' = * ) , T41 , 7 C T = '),T50,7( * = ' ) , T58 , 6 C ' = ') /) 

5030 FORMAT C1X,T6,F7.3,T15,F6 . 4, T23 , F7 .4,T32,E8.2,T41,F7.4,T50,F7 .4, 
+T58 , E9 . 3 ) 

5040 FORMAT C // 3 C IX, T5 , A17 , F8 . 4 /)) 

6 00 0 FORMAT C1X,F4.1,F1Q.3,F7.3) 

END 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

FUNCTION XFNCPHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

x THIS FUNCTION SUBPROGRAM IS USED TO DETERMINE STREAMWISE LENGTH 

x AS A FUNCTION OF THETA. XFN REPRESENTS PARAMETRIC RADIUS. 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL XFN, PHI , A, K, DEFF 
COMMON /GEOM/A , K, DEFF 

XFN = SQRTCCKXX2+.25X(CKXX2-1 .)XSINC2XPHI))XX2)/C (COSC PHI ) )xx2 
++CKXSINCPHI))XX2)) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION GC PHI ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

x THIS FUNCTION SUBPROGRAM REPRESENTS THE GRAVITY FUNCTION. 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL G, PHI, A, K, DEFF 
COMMON /GEOM/A, K, DEFF 

G = SINCPHI)/SQRT((SINCPHI))XX2+C KxCOS C PHI)) xx 2) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
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FILE: MIXED1 FORTRAN A 



FUNCTION Gl(PHI) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX3 exxxxxxxxxxxx 

X 

X THIS FUNCTION SUBPROGRAM REPRESENTS THE STREAMWISE DERIVATIVE 
x OF THE GRAVITY FUNCTION. 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
REAL G1,PHI,A,K,DEFF 
COMMON/GEOM/A , K, DEFF 

Gl=Kxx2xC0S(PHI)/SQRTCC(SIN(PHI))xx2+(Kxe0SCPHI))xx2)xx3) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION FI C PHI ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

* THIS FUNCTION SUBPROGRAM REPRESENTS THE VAPOR VELOCITY 

x FUNCTION 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL FI, PHI, A, K, DEFF 
COMMON /GEOM/A, K , DEFF 

FI = ( 1 .+K)XSINCPHI)/SQRT(CSINCPHI))XX2+CKXC0SCPHI))XX2) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION Fll(PHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X THIS FUNCTION SUBPROGRAM REPRESENTS THE STREAMWISE DERIVATIVE 
X OF THE VAPOR VELOCITY FUNCTION, 
x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL FIX, PHI, A, K, DEFF 
COMMON/ GEOM/A, K, DEFF 

Fll=(l .+K)XKXX2XC0S(PHI)/SQRT( (CSIN(PHI) )XX2+(KXC0S(PHI) )XX2)XX3) 
END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION F2CPHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X THIS FUNCTION SUBPROGRAM REPRESENTS THE PRESSURE GRADIENT 

x FUNCTION, 

x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL F2,PHI,A,K,DEFF 
EXTERNAL XFN 
COMMON/GEOM/A, K, DEFF 

F2=C (1 .+K)XJOXX2XSIN(2XPHI)/((SIN(PHI))XX2+(KXCOS(PHI))XX2) 

+/XFNC PHI ) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION F3CPHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X THIS FUNCTION SUBPROGRAM REPRESENTS THE SURFACE FUNCTION. 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
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FILE: MIXED I FORTRAN A 



REAL F3,PHI,A,K,DEFF 
EXTERNAL XFN 
COMMON/GEOM/A, K, DEFF 

F3 = K*(1 .-KXX2)XSIN(2XPHI)/SQRT(((SIN(PHI))XX2+(KXC0S(PHI))XX2) 
+xx5 )/XFN( PHI ) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION FCN(PHI,Z) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X THIS FUNCTION SUBPROGRAM IS USED IN A PROBLEM-SOLVER 
X SUBROUTINE TO SOLVE THE O.D.E. 
x DZ/DX = F( PHI / Z) 

x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL F21,F31,PHI,FCN1,FCN2,FCN,A,K,DEFF,F,B0,P 

COMMON /GEOM/A,K,DEFF/PARAM/F,BO,P 

EXTERNAL XFN , G, GI , FI , FU, F2, F3 

CALL DRV C F2, PHI , F21 ) 

CALL DRVCF3,PHI,F31) 

FCN1=C2.XFXG(PHI)+DEFF/AXPXF2CPHI)+3.X(DEFF/A)XX2XFXB0XF3(PHI)) 

+XZXX3+F1(PHI)XZ 

FCN2=2.XAXXFNCPHI)/DEFF-(2.XFXG1CPHI)+DEFF/AXPXF21+3.X(DEFF/A)XX2 

+xFxB0xF31)xZxxA/3.-F11CPHI)XZXX2 

FCN=FCN2/FCN1 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SUBROUTINE DRV( FN, PHI , FN1) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X THIS SUBROUTINE NUMERICALLY DIFFERENTIATES A FUNCTION. 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
REAL FN,PHI,FN1,EPS 

FN1=(FNCPHI+1 .E-6)-FN(PHI-l.E-6))/2.E-6 

RETURN 

END 
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APPENDIX C - COMPUTER CODE FOR TWO-PHASE BOUNDARY LAYER 



MODEL 

^ PROGRAM VAPBL1 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

* THIS PROGRAM NUMERICALLY SOLVES THE PROBLEM OF MIXED 

H CONVECTION CONDENSATION WITH SURFACE TENSION ON AN ELLIPTICAL 
x CYLINDER. THE PROBLEM FORMULATION USES THE SOLUTION OF THE 
X TWO-PHASE BOUNDARY LAYER EQUATIONS FOR THE CONDENSATE AND VAPOR. 

* 

x THE PROGRAM USES AN IMSL NUMERICAL PROBLEM SOLVER, DIVPAG, 

X TO SOLVE A 1ST ORDER ORDINARY DIFFERENTIAL EQUATION. THE 
x PROBLEM SOLVER IS WRITTEN IN DOUBLE PRECISION WHILE THE MAIN 
x PROGRAM IS IN SINGLE PRECISION. THE MAIN PROGRAM MUST BE 

* COMPILED IN DOUBLE PRECISION CAUTO DOUBLE). 

X 

* THIS PROGRAM CAN BE RUN WITH DIMENSIONLESS PARAMETERS 

X (F,G) SPECIFIED OR WHERE THE PARAMETERS ARE CALCULATED 

x BASED ON SPECIFIED FLUID PROPERTIES. MODIFICATIONS TO THE 

x PROGRAM (REMOVAL OF UNNECCESSARY STEPS BY "COMMENTING" THE 

x STATEMENT OUT OF THE PROGRAM) IS REQUIRED. SPECIFIED FLUID 

x PROPERTIES INCLUDE T SAT, T WALL, AND UINF. 
x 

X IF THE ALGORITHM HAS DIFFICULTY LOCATING THE VAPOR SEPARATION 

* POINT, THEN THE PROGRAM NEEDS TO BE MODIFIED TO INPUT THE 

* PARAMETRIC AND, PHI, FOR THE SEPARATION POINT, 
x 

X NUMERICAL RESULTS ARE DIRECTED TO AN EXTERNAL DATA FILE. 

X 

x MAJOR VARIABLES USED ARE: 
x 



X 


PHI 


PARAMETRIC ANGLE 


X 


X 


STREAMWISE LENGTH 


X 


DEFF 


EFFECTIVE DIAMETER 


X 


DEL 


CONDENSATE FILM THICKNESS 


X 


DELND 


DIMENSIONLESS FILM THICKNESS 


X 


2 


DIMENSIONLESS VAPOR MOMENTUM THICKNESS 


X 


F, JARPRL 


DIMENSIONLESS PARAMETERS 


X 


UINF 


VAPOR FREE STREAM VELOCITY 


X 


TS 


VAPOR SATURATION TEMPERATURE 


X 


TW 


WALL SURFACE TEMPERATURE 


X 


DENSL 


CONDENSATE DENSITY 


X 


DENSV 


VAPOR DENSITY 


X 


CONDL 


CONDENSATE THERMAL CONDUCTIVITY 


X 


VISL 


CONDENSATE VISCOSITY 


X 


VISV 


VISCOSITY OF VAPOR 


X 


HFG 


LATENT HEAT OF VAPORIZATION 


X 


NU 


NUSSELT NUMBER 


X 


REV 


VAPOR REYNOLDS NUMBER 


X 


RE 


TWO-PHASE REYNOLDS NUMBER 


X 


NUAVG 


MEAN NU NUMBER 


X 


ALPHA 


HEAT TRANSFER COEFFICIENT 


X 


QFLUX 


HEAT FLUX 


X 


A 


SEMI-MAJOR AXIS 


X 


K 


ECCENTRICITY 


X 


XFN 


PARAMETRIC RADIUS 


X 


G 


GRAVITY FUNCTION 


X 


FI 


VELOCITY FUNCTION 


X 


Fll 


DERIVATIVE, VELOCITY FUNCTION 


X 


F2 


PRESSURE GRADIENT FUNCTION 


X 


F3 


SURFACE TENSION FUNCTION 


X 


TAU 


INTERFACIAL SHEAR 


X 


KAPPA 


PRESSURE GRADIENT FORM PARAMETER 


X 


KAPPAA 


SUCTION FORM PARAMETER 



x 



xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

PARAMETERC PI =3 . 141592654, GRA V = 9. 807, NEQ = 2, NEQ 1=1, NP A RAM = 50) 

REAL PHI( 0: 180000 ),DEL(0: 180000 ),X<0: 180000) ,A,TS,TW, 
+H,A1,A2,A3,A < 4,A5,A6,A7,A8,A9,A10,A11,TF, PSAT ,TB,TC,TD,R,TE,DENSV, 
+TL, DENSL, VI SL, VISV, TG, CONDL , HFG, SIGML , F,K,X1 ,X2,X3, DEFF,UINF, 

+ REV , RE, JARPRL , DELOO , DELO, ZO , TAU, DTAU, FDEL 0 , DFDELO , DELI , Z1 , ALPHAO , 
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oo oo 



FILE: VAPBL 



FORTRAN A 



+SUM, AAC 1 , 1 ) , FCN, FCNJ , HI NIT, PA RAM ( NPARAM) , KAPPA, KAPPAA , DKAPPA, 
+T0L,XEND,X3EG, Y(2),DDEL, DZ, Y1 ( 1 ) , NU , NUAVG, NUORE, DKAPPAA , XO , DXO, 
+DELND, DELL , ALPHA, QFLUX, DQFLUX, D2UPHI, DUPHI , UPHI 

INTEGER N, J J , I I , IDO , IMETH 

COMMON /GEQM/A , K, DEFF/PARAM/ J ARPRL , F 

EXTERNAL ZFN , DELFN , XFN, G, G1 , FI , FI 1 , NUSSELT , DRV , FCN, FCNJ , IVPAG, 
+SSET 

0PEN(20,FILE=»/VAP0R1 OUTPUT A1 ' , STATUS 3 f OLD* ) 

OPEN( 30 , FILE= '/VAP0R2 OUTPUT A1 ’ , STATUS= 'OL D’ ) 

PRINT*, * ELLIPSE OF THE FORM (X/B)**2+(Y/A)**2=1 ' 

PRINT*, ' 1NPUT "A" DIMENSION CM)? 1 
READ*, A 

PRINT*, 'INPUT T SAT AND T WALL 00?’ 

READ*, TS, TW 

PRINT*, 'INPUT ANGLE STEP SIZE (DEGREES)? 1 
READ*, H 

PRINT*, 'STEP SIZE, H = ' , H 
PRINT*, 'INPUT PHI AT SEPARATION?' 

READ*, PHIS 

PRINT*, 'INPUT UINF?' 

READ*, UINF 
PRINT*, 'INPUT K?* 

READ*, K 

PRINT*, 'INPUT G?' 

READ*, JARPRL 

N=180/H+l 
PRINT*, 'N-1=',N-1 
H=PI/180.*H 

*** DETERMINE FLUID PROPERTIES ******** 

A1 =15 . 49217901 
A2=-5. 6783717693 
A3 = 1 . 45975846 37 
A4=13. 877000608 
A5=-80. 887673591 
A6=123. 56883468 
A7=-188 .321212064 
A8=660. 91763485 
A9=-1382. 4740091 
A10=1300. 1040184 
All=-449 .39571976 
TF=TS/ 1000. 

PSAT = 1 . E+6*EXP( A1+A2/TF+A3*L 0G(TF)+TF*A4+A5*TF**2+A6*TF**3+A7* 
+TF**4+A8*TF**5+A9*TF**6+A1 0*TF**7 +A1 1*TF**8 ) 

TB=1500./TS 

TC=2 . 5*L0G( 1 . “EXP( -TB) ) 

TD=. 001 5/(1.+. 0001 *TS)~. 000942*EXP ( TB+TC) *SQRT (l./TB)-.0004882*TB 
R=46l .51 

TE=2.*PSAT/(R*TS) 

DENSV=TE/ ( 1 . +SQRT ( 1 . +2 . *TD*TE) ) 

TL=(2.*TW+TS)/3. 

DENSL=l./( .0012674 -TL*(2. Q2915E-6-TL*3 . 8333E-9 ) ) 
VISL=2.414E-5*10**(247.8/(TL-140. )) 

VISV=-4.478415E-6+TS*(5. 0216E-8-1 .579E-11*TS) 

TG=TL/27 3.15 

C0NDL=-.92407+TG*(2.8395-TG*(l .8007-TG*( . 52577-TG* . 07344) ) ) 
HFG=3468920.-TS*(5707 .4-TS*(ll .5562-TS*. 0133103)) 
SIGML=(-.0003*(TL-273.15)**2-.138*(TL-273.15)+75.6)/1000. 

*** DETERMINE PHI ( I ) , X( I ) , DEFF ******** 

*** USES SIMPSON'S RULE TO NUMERICALLY INTEGRATE ALONG SURFACE *** 

PHI ( 0 ) =0 . 

X(0)=0. 
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FILE: VAPBL 



FORTRAN A 



DO 10, 1=1, N 

PHICI)=PHICI-1)+H 

X1=XFNCPHICI-1)) 

X2 S XFNC C PHI C 1-1 )+PHI C I ) )x . 5) 

X3 = XFN( PHI CD) 

XC I ) =X( 1-1 )+H/6 . ^CXl+4 . XX2+X3) 

10 CONTINUE 

DEFF=2.XXCN)/PI 

PRINTx, 'COMPLETED DE DETERMINATION, DE= 1 , DEFF 
xxxx DIMENSIONLESS PARAMETER DETERMINATION xxxx 

F=VISL*DEFF*HFG/(C0NDLxUINFXX2x(TS-TW)) 

REV=UINFXDENSVXDEFF/VISV 

RE=UINFXDENSLXDEFF/VISL 

JARPRL=C0NDLXCTS-TW)/HFGXSQRTC DENS L/C DENSVXVISL XVI SV) ) 

PRINT*, COMPLETED PARAMETER DETERMINATION 1 

WRITEC 20 ,50 00) 'SAT. TEMP. CK) = ',TS,'WALL TEMP. = 1 , TN, 

+ 1 F = 1 , F, 1 UINF CM/S) = 1 , UINF, 'G = ' , JARPRL 
WRITEC 20,5010) 1 "A" LENGTH CM) = 1 , A, 1 "K" ECCENTRICITY = »,K, 

+ T DEFF CM) = ' , DEFF 

xxx DETERMINE INITIAL VALUES OF Z AND DEL **X 

xxx USES INTIAL VALUE OF DEL FROM SHEKRI L ADZE-GOMEL AURI xxx 

xxx FOR MIXED CONVECTION AS STARTING POINT TO FIND Z THEN xxx 

xxx REDETERMINES DEL FOR THE PREDICTED VALUE OF Z USING xxx 

xxx NEWTON 1 S METHOD, CHECKS FOR CONVERGENCE AND ITERATES xxx 

DEL00=SQRTC . 75*K/Fx( SQRTC C C 1 . +K)/K) XX2+16 . xF*A/ C 3 . XDEFF) ) 

+-C 1 . +K)/K) ) 

DEL0=DEL00 

Z0=C-. 0735XXFNC0. )XJARPRL/ DEFF/ F11C 0. )/DELO+SQRTCC .0735* 

+XFNCO . )XJARPRL/DEFF/F11C0. )/DEL0)xx2+. 0735XXFNC0. )/DEFF/ 

+ F11 CO.)) )xx2 

UPHI = F1 C 0 . ) 

DUPHI =F11 CO.) 

CALL DRVCF11, 0 . ,D2UPHI) 

X0=XFN CO.) 

CALL DRVCXFN, 0 . , DXO) 

20 KAPPA=DEFF/XQXDUPHIXZO 

30 KAPPAA=.0682+.174XJARPRL*SQRTCZ0)/DEL0 

DKAPPA=DEFFxZ0xD2UPHI/X0-DEFFxDUPHIxZ0xDX0/X0xx2 
DTAU=6 .44XSQRTCKAPPAXKAPPAA+KAPPAAXX2)XDUPHI/SQRTCZ0)+3.22XUPHI 
+XKAPPAAXDKAPPA/SQRTCZQ)/SQRTCKAPPAXKAPPAA+KAPPAAXX2) 

FDEL0=X0/DEFF-F/3.XG1C0. )XDEL0XX4-DTAUXDEL0XX3/C4.XJARPRL) 
DFDEL0=-4. XFXG1 C 0 . )XDELQxx3/3 . -3 . XDTAUXDEL0XX2/C 4 . x JARPRL ) 

DELI = DEL 0-FDEL 0/ DFDEL 0 

IFC C ABSC DEL 1-DEL 0) . GT . C . 001XDEL 0 ) ) . AND . C FDEL 0 . GT . 1 . E-7 ) ) THEN 
DEL Q = DEL1 
GO TO 30 
ENDIF 

Z1=C-. 0735XXFNC0. )xj ARP RL/ C DEFFXF1 1 C 0 .) XDEL 1 )+SQRTC C . 0735* 

+XFNC0 . )XJARPRL/CDEFFXF11C0. )XDEL1))XX2+.0735XXFNC0. )/CDEFFX 

+F11C0. ))))XX2 

IF CCABSCDEL1-DEL00) .GT. C .OOIXDELOO)) . OR. CABSCZ1-Z0) .GT. C .001* 
+Z0))) THEN 
Z0 = Z1 

DEL 00 = DEL 1 
DEL0=DEL1 
GO TO 20 
ENDIF 

DELC 0 ) =DEL1 
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FILE: VAPBl 



FORTRAN A 



xxx 

XXX 



c 



XXX 

100 

190 



XXX 

XXX 

XXX 



200 



XXX 



ALP HA 0=1. /D EL (0) 

5UM=Q . 

PRINTS, ’COMPLETED ZC0),DELC0) DETERMINATION 1 ,* Z=*,Z1, 

+ * DEL = * , DEL CO) 

EVALUATION OF ZCI) AND DEL ( I) **** 

USES IMS L ROUTINE DIVPAG **** 

HINIT = l.E-<+ 

IMETH=2 

CALL SSETCNPARAM,Q.0,PARAM,1) 

PARAMC 1 ) =HINIT 
PARAMC12)=IMETH 
XBEG=PHI CO) 

Y( 1 ) =Z1 
Y( 2) = DEL CO) 

T0L=l.E-3 
IDO = l 

DO 100 , 1=1 , N-l 

IFCPHICI) .GT.PHIS) GO TO 190 
XEND=PHI ( I ) 

CALL D I VP AG ( IDO , NEQ, FCN , FCNJ , AA , XBEG, XEND, TOL , PARAM, Y) 

DEL C I ) =YC 2 ) 

KAPPA=DEFF/XFNCPHICI))*F11CPHI(I))*Y(1) 

KAPPAA=.0682+.174*JARPRL*SQRTCYC1))/YC2) 

TAU=6 . 44*F1CPHICI))*SQRTCKAPPA*KAPPAA+KAPPAA**2)/SQRTCYC1)) 
WRITEC 30 , X) »PHI=*,PHICI)*180/PI, *TAU=* ,TAU 
II = I 



CHECK FOR VAPOR BOUNDARY LAYER SEPARATION *XXX 

IFCTAU.LE.O.) GO TO 190 
CONTINUE 

PRINT*, ’SEPARATION ANGLE AT PHI= ' , PHI C II )X180 ./PI 
PRINT*, *TAU= * , TAU, * DEL= * , DEL C I I ) 

ID0 = 3 

CALL DIVPAGC IDO, NEQ, FCN, FCNJ, A A, XBEG, XEND, TOL , PARAM, Y) 

USES A SIMPLE NUSSELT ANALYSIS TO COMPUTE FILM THICKNESS X*** 
DOWNSTREAM OF THE SEPARATION POINT. PROBLEM SOLVER IS x*** 

IMS L DIVPAG **** 

IFCII.NE.CN-1)) THEN 
XBEG=PHI C II ) 

IDO = l 

Y1 C 1) =DEL (II) 

DO 200, 1=11+1, N-l 
XEND=PHI ( I ) 

CALL DIVPAGC IDO, NEQ1, NUSSELT, FCN J, AA, XBEG, XEND, TOL, PARAM 
+ ,Y1) 

DEL Cl) =Y1 Cl) 

CONTINUE 
ENDIF 
ID0 = 3 

CALL DIVPAGC IDO, NEQ1, NUSSELT, FCN J,AA, XBEG, XEND, TOL, PARAM, Yl) 

PRINT*, *COMPLETED ZCI), DELCI) DETERMINATION* 

WRITEC 20 , 5050 ) PHI C II )*180 ./PI 
WRITEC 20 , 5020 ) 

DETERMINE HEAT TRANSFER PROPERTIES ****xxxxx 



JJ = 0 

DO 300, 1=1, N-l 
DELND = DEL Cl) 

DEL L = DELND*SQRTC VISL*DEFF/C DENSL*UINF) ) 

ALPHA = 1 ./DEL Cl) 

SUM = SUM+C ALPHA 0*XFNC PHI CI-1))+ALPHA*XFNCPHICI)))*H*.5 
ALPHAO=ALPHA 
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FILE: VAPBL 



FORTRAN A 



IFCCJJ.EQ.I) .OR.CJJ.EQ.O)) THEN 
QFLUX=ALPHA*CTS-TW) 

NU=DEFF/DELL 
DQFLUX- 1 ./DELND 
JJ=JJ+10./130.*N+1 

WRITEC 20 , 5030 ) PHI ( I )X180/PI , X( I)*2/ ( DEFF*PI) , DELL 
+ *1000 ,QFLUX/ 1 000, DELND, DQFLUX, NU 

ENDIF 

300 CONTINUE 

WRITEC 20 , 5030) PHI ( N-l ) *130 . /PI , XC N-l ) *2./ ( DEFF*PI ) , DEL L*1 00 0 . , 
+QFLUX/1 000 . , DEL ND, DQFLUX, DEFF/DELL 
WRITEC 30,*) XCN-1)*2./CDEFF*PI),DELND 

NUAVG=SUM*2/CPI*DEFF) 

WRITEC2Q , 5040 ) ' NU/SQRT C RE) = *, NUAVG/SQRTC RE) 

5000 FORMAT C1X,T5,7QC '* ' ) // IX, T5, 'HEAT TRANSFER PROPERTIES' // 2C1X, 
+T1 0 , A17 , F6 . 2 / ) ,5C1X,T10,A17,E9.3 /)) 

5010 FORMAT C / IX, T5, ' ELLIPSE GEOMETRY' // 3C IX, T1 0 , A19 , F6 . A / )) 

5020 FORMAT C / IX, T7 , 'THETA ', T15, ' DIM X ' , T25 , ' DEL ' , T32 , ' HT FLUX', 

+TA1 , ' DIM DEL',T52, 'DIM',T58, 'NUSS #' / T7 , ' C DEG) » , T25, ' CMM) ' , T32, 
+ » C KW/M2) ' , T50 , 1 HT FLUX' / T6 , 7 C ' = ' ) , T15, 6 C ' = » ) , T23, 7 C ' = ' ) , T32, 

+7C »='),TA1,7C '='),T5Q,7C'='),T53,6C ' = ») /) 

5030 FORMAT C1X,T6,F7.3,T15,F6.A,T23,F7.A,T32,E8.2,TA1,F7.A,T50,F7.A, 
+T58 , F9 . 1 ) 

50 AO FORMAT C // 3 C IX, T5 , A17 , FS . A /)) 

5050 F0RMATC1X,T6, 'VAPOR SEPARATION AT PHI = »,F6.2) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION XFNCPHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

* THIS FUNCTION SUBPROGRAM IS USED TO DETERMINE STREAMWISE LENGTH 

* AS A FUNCTION OF PHI. 

x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL XFN, PHI , A, K, DEFF 
COMMON /GEOM/A, K, DEFF 

XFN=A*SQRTCCK**2+.25x(CK**2-l.)*SINC2*PHI))**2)/CCC0SCPHI))**2 
++C KxSINCPHI ) )**2) ) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION GCPHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

* GRAVITY FUNCTION 
X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL G, PHI , A, K, DEFF 
COMMON /GEOM/A, K, DEFF 

G=SINCPHI)/SQRTCCSINCPHI)) **2+C K*COS C PHI ) ) **2 ) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION G1CPHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

* STREAMWISE DERIVATIVE OF GRAVITY FUNCTION 
X 
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FILE: VAP3L 



FORTRAN A 



XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

REAL G1,PHI,A,K,DEFF 
COMMON/GEQM/A/K/DEFF 

G1=KXX2XC0S(PHI)/SQRT(((SIN(PHI))XX2+(KXC0S(PHI))XX2)XX3) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION FI (PHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

x VAPOR VELOCITY FUNCTION 
X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL FI / PHI / A, K, DEFF 
COMMON /GEOM/A,K,DEFF 

FI —Cl . +K)XSIN(PHI)/SQRT( (SIN(PHI))XX2+(KXC0S(PHI))XX2) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION FI 1 (PHI ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X STREAMWISE DERIVATIVE OF VELOCITY FUNCTION 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL F11,PHI,A,K,DEFF 
COMMON/GEOM/A, K, DEFF 

FI 1 = ( 1 .+K)XKXX2XCQS(PHI)/SQRT(((SIN(PHI))XX2+(KXC0S(PHI))XX2)XX3) 
END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION DELFN( PHI , Z, DEL / DZ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X STREAMWISE DERIVATIVE OF CONDENSATE FILM THICKNESS 
X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL DELFN, PHI , Z, DEL , DZ, DELFN1 , DELFN2, DELFN3, X, FF1 / FF1 1 / FF12, DXFN 
+/A/K, DEFF, JARPRL/F/ KAPPA, KAPPAA/ XX/ FKAPPA 

EXTERNAL XFN/ FI / Fll /G/Gl 



COMMON/GEOM/Az K/ DEFF/PARAM/ J ARPRL / F 

IF( PHI . EQ . 0 . D+O ) THEN 
DEL FN = 0 . D+O 
ELSE 

FF1 =F1 ( PHI ) 

FF1 1 S F1 1 ( PHI ) 

CALL DRV( Fll /PHI, FF12) 

X=XFN( PHI ) 

CALL DRV(XFN,PHI/ DXFN) 

KAPPA=DEFFXFF11XZ/X 
KAPPAA=.0682+.174XJARPRLXSQRT(Z)/DEL 
XX=KAPPAAXKAPPA+KAPPAAXX2 
FKAPPA=SQRT (XX) 

DELFNl=.S05XKAPPAAXDEFFXFFlxFFllXDELXX3/( JARPRLXFKAPPAXXXSQRT(Z)) 
++(7 . 0035E-2)X(2.XKAPPAA+KAPPA)XFF1XDELXX2/(FKAPPAXZ)-.S05XFKAPPAX 
+FFlXDELXX3/( JARPRLXZXX1 .5) 
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FILE: VAPBL 



FORTRAN A 



DELFN2=F*GCPHI)*DEL*x3+3. 22*FF1*FKAPPA*DEL **2/ ( JARPRLXSQRTC Z) ) 
+-.14007X(2.XKAPPAA+KAPPA)XFF1XDEL/FKAPPA 

DELFN3=X/DEFF-FXGl(PHI)XDELXX4/3.-1.6lXFKAPPAxFFllXDELXX3/(JARPRL 
+XSQRTCZ) )- . 305XKAPPAAXDEFFXFF1XFF12XSQRT ( Z) XDELXX3/ ( J ARPRLxFKAPPAx 
+X)+.305XKAPPAAXDEFFXDXFNXFF1XFF11XSQRT(Z)XDELXX3/C J ARPRLXFKAPPAX 
+XXX2) 

DELFN=CDELFN3-DELFN1XDZ)/DELFN2 

END IF 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION ZFNC PHI ,Z, DEL ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X STREAMWISE DERIVATIVE OF VAPOR MOMENTUM THICKNESS PARAMETER 
X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL ZFN,PHI,Z,DEL,DEFF,A,K, JARPRL,F 

EXTERNAL XFN 
EXTERNAL FI 
EXTERNAL Fll 

COMMON/GEOM/A,K,DEFF/PARAM/JARPRL,F 

IFCPHI.EQ.O.D+O) THEN 
ZFN=0 . D+O 
ELSE 

ZFN=( ,44lx(l .-2 .xJARPRLxSQRT(Z)/DEL)-6,xDEFFXF11CPHI)xZ/XFN(PHI)) 
+XXFNC PHD/F1CPHD/DEFF 
ENDIF 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE NUSSELTCNEQ, PHI , Y, YPRIME) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X STREAMWISE DERIVATIVE OF CONDENSATE FILM THICKNESS USING 
X A NUSSELT ANALYSIS DOWNSTREAM OF VAPOR SEPARATION POINT, 
x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INTEGER NEQ 

REAL PHI,YCNEQ),YPRIME(NEQ),A,K,DEFF, JARPRL,F 
EXTERNAL XFN,G,G1 

COMMON/GEOM/A, K, DEFF/PARAM/ J ARPRL , F 

YPRIMEC 1 ) =(XFN(PHI )/DEFF-F/3 . XYC 1 )XX$xGl (PHI ) )/( FXGCPHI )xY( 1 ) xx3) 

RETURN 

END 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX* 

SUBROUTINE DRVC FN, PHI , FN1 ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

x SUBROUTINE TO NUMERICALLY DIFFERENTIATE A FUNCTION 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
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FILE: VAPBl 



FORTRAN A 



REAL FN , PHI , FN1 , EPS 

FN1 = C FN( PHI + 1 . E-6 ) -FNC PHI-1 . E-6 ) )/2 . E-6 

RETURN 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SUBROUTINE TAUFNC PHI , Z, DEL , DZ, DDEL , TAU , DTAU ) 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

x SUBROUTINE TO DETERMINE THE INTERFACIAL SHEAR STRESS AND 
x ITS STREAMWISE DERIVATIVE. 

X 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

REAL PHI , Z, DEL , DXFN, DZ, DDEL , KAPPA , KAPPA1 , KAPPA A, DKAPPA, 

+DKAPPA1 , DKAPPAA, TAU, DTAU, A , K, DEFF, JARPRL , F 

EXTERNAL XFN,F1,F11 

COMMON/GEOM/A, K, DEFF/PARAM/ JARPRL , F 

CALL DRV ( FI 1 , PHI , F12) 

CALL DRV (XFN, PHI , DXFN) 

KAPPA=DEFFXF11(PHI) XZ/XFNC PHI ) 

KAPPA1=JARPRLXSQRT(Z)/DEL 

KAPPAA=.0682+.174*KAPPA1 

DKAPPA=DEFF/XFNC PHI )x (ZXF12+F11 ( PHI) XDZ) -DEFF/ CXFNC PHI ))XX2X 
+FII(PHI)*ZXDXFN 

DKAPPA1= JARPRL x( . 5XDELXDZ-ZXDDEL )/ ( SQRTCZ) XDELXX2) 
DKAPPAA=.17AXDKAPPA1 

TAU = 6 .44xFIl(PHn*SQRT((KAPPAAxKAPPA+KAPPAAxX2)/Z) 

DTAU=6 . 44 x( C (KAPPAAXKAPPA+KAPPAAXX2) XF1 1 ( PHI ) + . 5xFl ( PHI )X 
+CKAPPAAXDKAPPA+C2.XKAPPAA+KAPPA)XDKAPPAA))/SQRT(ZXCKAPPAAX 

+KAPPA+KAPPAAXX2))-.5X(KAPPAAXKAPPA+KAPPAAXX2)XF1(PHI)XDZ/ 

+SQRTC ZXX3XC KAPPA AXKAPPA+KAP PA Axx2) ) ) 

RETURN 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SUBROUTINE FCNCNEQ, X, Y, YPRIME) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

x SUBROUTINE TO INPUT DIFFERENTIAL EQUATION TO BE SOLVED 
X BY THE PROBLEM SOLVER. 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INTEGER NEQ 

REAL X,Y(NEQ) , YPRIMEC NEQ) 

EXTERNAL ZFN, DELFN 

YPRIMEC 1 ) =ZFN(X, Y( 1 ) , YC 2) ) 

YPRIMEC 2) = DELFN(X/ YC 1) , YC 2) /YPRIMEC 1) ) 

RETURN 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SUBROUTINE FCNJ C NEQ, X, Y, DYPDY ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
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>K W * 



FILE: VAPBL 



FORTRAN A 



SUBROUTINE REQUIRED BY PROBLEM SOLVER, NO FUNCTION IN ANALYSIS 
INTEGER NEQ 

REAL X,Y(NEQ),DYPDYOO 

RETURN 

END 
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